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WORK STATEMENT 


The co-operative agreement was initiated for achieving the following: 


1. Carry out a tip-to-tail computation of the flow field of a generic hypersonic vehicle in the 
power-on configuration. The main focus of the task was to establish and demonstrate a solution 
procedure/strategy that would be useful for solving the entire flow field associated with an 
airbreathing hypersonic vehicle. The plan was to utilize available Computational Fluid 
Dynamics (CFD) codes for this purpose rather than generating a new specialized CFD solver. 
This would need modifications to the codes to enhance their capabilities in terms of accurate 
models for physical features such as turbulence, finite rate chemistry, improved boundary 
conditions etc. 


2. Develop and demonstrate a turbulence-chemistry interactions model applicable for high- 
speed reacting flows involving finite rate chemistry. 

The following sections describe the milestones achieved, work in progress and 
suggestions for future work pertaining to the items described above. In this report, 
bibliographical references to the physical models such as the turbulence models, thermodynamic 
models, finite rate chemistry models etc. and CFD codes used are omitted. However, these 
references can easily be obtained from the attached publications pertaining to the present task. 


HYPERSONIC TIP-TO-TAIL COMPUTATION 


The geometry of the generic hypersonic vehicle configuration was provided by Dr. Tom 
Edwards of the Computational Aerosciences Branch, Fluid Dynamics Division at NASA Ames 
Research Center. The geometric data included body definitions and preliminary grids for the 
forebody (nose cone excluded), midsection (propulsion system excluded) and afterbody sections. 
This data was to be augmented by the nose section geometry (blunt conical section mated with 
the non circular cross section of the forebody initial plane) along with a grid and a detailed 
supersonic combustion ramjet (scramjet) geometry (inlet and combustor) which should be 
merged with the nozzle portion of the afterbody geometry. The solutions were to be obtained by 
using a Navier-Stokes (NS) code such as TUFF for the nose portion, a Parabolized Navier 
Stokes (PNS) solver such as the UPS and STUFF codes for the forebody, a NS solver with 
finite rate hydrogen-air chemistry capability such as TUFF and SPARK for the scramjet and a 
suitable solver (NS or PNS) for the afterbody and external nozzle flows. Proper interfacing was 
to be maintained between the solutions of the different sections. 

NOSE SECTION 

The nose is assumed to be a blunt one in the present computations in keeping with the 
accepted designs for an airbreathing hypersonic vehicle. As mentioned above, the nose section is 
a blunt-nosed cone at zero angle of attack. An axisymmetric grid was generated for the nose 
section using the graphic utility VG available on workstation win35. The TUFF code was used 
to solve the flow field upto the axial location corresponding to the first plane of the forebody 
geometry. Finite rate nonequilibrium air-chemistry model was used along with a two-equation 
turbulence (k-epsilon) model. The flow Mach number was 15 and the Reynolds number based on 
unit length was 1,700,000. The surface of the body was assumed to be maintained at a constant 
temperature of 243.4 K. No-slip boundary conditions were used at the body surface. The 
solutions were then interpolated onto the first plane of the forebody grid. 

FOREBODY 

Computations were carried out using the PNS solvers STUFF and UPS in order to 
ascertain that these codes were equivalent as for as the present application was concerned. The 
nose-section solution was used along with the forebody grid for this purpose. The solutions were 
obtained for the flow field around the entire forebody upto the inlet location. Baldwin-Lomax 
turbulence model was used since the version of the UPS code used does not have the two- 



equation turbulence model capability. The solutions were compared (temperature, velocities and 
pressure) and found to be in excellent agreement. Since the TUFF and STUFF codes have 
similar structure and physical models, it was decided to use them for the NS and PNS 
applications, respectively, unless conditions dictated the utilization of a different code such as 
SPARK. This would be extremely useful in interfacing solutions between different sections of 
the vehicle. 

Since the nose section was a conical body with circular cross section, it was necessary to 
generate a series of cross sections at designated intervals which changed in shape from a circle to 
that of the first cross section of the forebody grid. A three dimensional grid that bridged the last 
plane of the nose section with the first plane of the forebody was then generated using these 
surface definitions (cross sections) and was combined with the given forebody grid to form the 
new three dimensional grid for the forebody flow field solutions. These solutions were obtained 
by using the space marching PNS solver STUFF. The nose-section solutions were used as initial 
data for this purpose. The two-equation (k-epsilon) turbulence model including low-Reynolds 
number terms to account for the near-wall flow field was used and the finite rate air-chemistry 
model of the nose section was maintained. Even though the turbulence field was not one of the 
primary targets in these computations, turbulence modeling at the two-equation level was 
maintained in anticipation of the inlet/combustor flow simulations which require the use of the 
physically more realistic two-equation model rather than the simple and easy to use algebraic 
Baldwin-Lomax model. As a result, the computations required more computer resources. 
Typical results (composite since the full computational data at all the x-locations could not be 
stored) are shown in the figures at the end of the report. Computationally, there were no 

unforeseen difficulties in obtaining these solutions. This task was completed in September, 
1992. 

INLET 

The inlet computations represent the critical (bottleneck) part of the entire tip-to-tail 
solutions. The difficulty was associated with the proper choice of a realistic inlet configuration 
at the same time minimizing the computational effort which utilized the available resources in an 
efficient manner. The geometric restrictions imposed by the limited length of the entire scramjet 
propulsion system to achieve the power requirements for the hypersonic flight is a major factor 
to be considered here. The mixing between the incoming airflow and the injected fuel dominates 
the flow field and the mixing efficiency plays a major role in determining the dimensions and 
geometry of the inlet. Also, physical aspects of the flow such as inlet unstart, flow spillage at 



the cowl location to maintain the proper air flow etc. are key factors that must be dealt with in the 
solutions. The solutions are computationally intensive, requiring an enormous amount of 
computational resources. Since the incoming forebody flow is three dimensional, the inlet flow 
solutions also require fully three dimensional analyses. 

The computational effort was divided into a series of tasks. First, an inlet 
configuration had be selected. Since the present work was meant to be a demonstration type with 
the purpose of establishing a solution procedure/strategy for carrying out tip-to-tail 
computations, it was decided to use a multi-module scramjet inlet for the present task. This 
would conform to the widely accepted inlet configurations being used in research and 
development efforts earned out at a majority of establishments working in the area of hypersonic 
flight. In order to reduce the computational effort thus minimizing the demand for computational 
resources, a three-module configuration that spans the entire undersurface of the vehicle was 
chosen for the present task. The next task was concerned with the selection of the geometric 
configuration of each inlet module. Once again, taking into account the realistic configurations 
and keeping with the aim of minimizing computational effort, a swept-back twin strutted module 
with side-wall compression was chosen for the present task. A sweep angle of 30 degrees was 
chosen along with a strut angle of 6 degrees. A schematic view of the inlet module is given in 
the figures. 

The above mentioned choice for the inlet geometry immediately introduced a new 
problem. This problem was associated with the choice of the particular CFD solver to be used 
for the solutions. The multiple-module configuration with the cowl plate located downstream of 
the strut tip introduced a variety of boundary conditions at different locations in the flow field. 
For example, a typical module has solid walls on three of the six boundary surfaces that define 
its geometry (omitting for the time being factors such as injection at solid walls). The inlet and 
exit boundaries pose no serious problems also. However, the fourth boundary that is aligned 
with the cowl plate surface involves both no-slip and flow-through conditions which are not easy 
to implement in the CFD codes under consideration. Also, due to the presence of the struts 
inside the flow region, an enormous grid (single, 3-D) would be required to have any chance at 
obtaining a realistic solution. The grid needed would be impractical to be used with a NS solver. 
As a result, a multi-block NS solver that includes all the physical models used for the computing 
the flow field around the forebody was needed. The multi-block version of the TUFF code 
became available at that time and hence was selected as the CFD solver for the inlet flow 
analysis. However, the problems associated with mixed boundary conditions, alluded to earlier, 
still remained requiring modification of the multi-block TUFF code to account for mixed 
boundary conditions. These modifications were incorporated in the solver which required 



considerable resources in terms of manpower as well as computational resources for validation. 
The changes led to a system by which different boundary conditions (from among the original 
options provided in the solver) could be specified at every grid point on the six boundary 
surfaces. 

The fuel of choice for the present case is hydrogen. However, the chemistry model in the 
TUFF code was for hydrocarbon-air combustion. This led to further modifications to the code 
whereby a hydrogen-air (finite rate) chemistry model involving nine species (hydrogen, oxygen, 
water, OH, H, O and nitrogen) and 20 reaction steps (suggested by the NASP Technical 
Committee) was incorporated. Also, thermodynamic models valid over a wider range of 
temperature than before was incorporated (by the author of the code). Even though chemical 
reactions involving combustion are not expected in the inlet, the chemistry model is included for 
the inlet flow solutions so as to be consistent with the combustor solutions which would follow. 

The solution domain is divided into adjoining blocks/zones with distinct boundary 
conditions. This is done in such a way that the occurrence of mixed boundary conditions along 
the boundaries is kept to a minimum. A multi-block grid involving 7 blocks was generated for 
the solutions. A schematic of the block arrangement is shown in the figures along with a 
representation of the grid system. The solution of the flow in the inlet represents the most 
computationally intensive part of the entire solution process. The flow is both external and 
internal involving strong and weak shocks/expansions in the domain. The precompressed air 
from the forebody (undersurface, mainly) enters the multiple module inlet and undergoes further 
changes due to the geometry of the inlet. There is flow spillage at the cowl requiring to pay 
extra attention to the grid there. 

Solution of the inlet flow was carried out over a long period of real time. One of the 
reasons for this was the fact that the computer resources ran out in 1992 about four months 
before the start of the next account year, virtually stopping the computationally intensive inlet 
flow solutions for about four months. Only peripheral work could be done at that time. The 
solutions were finally completed by July 1993. Representative results are included at the end of 
the report. 

COMBUSTOR 

The outflow from the inlet were to be used for the solution of the flow field beyond the 
cowl lip location (all internal flow, wall bounded ) into the combustor. The grid for this section 



(throat+combustor) is already generated. In order to save computational resources, it was 
decided to carry out the solution corresponding to only one of the modules (rather than all the 
three). However, the grid for all the three modules has already been generated in case a full 
configuration solution is needed. The throat flow solution is completed and is to be used for the 
combustor flow solution. 

The combustor geometry involves a backward facing step at the end of the strut wall 
extension in order to facilitate injection of hydrogen into the air flow. A 7-species, 7-reaction 
steps finite rate chemistry model will be used for the combustor flow simulations. The 
turbulence model is still the two-equation model described before. 


TURBULENCE -CHEMISTRY INTERACTION MODEL 

The numerical simulation of the hypersonic propulsion system for the generic hypersonic 
vehicle is the major focus of this entire work. Supersonic combustion ramjet (scramjet) is such a 
propulsion system as mentioned in the previous section. Hence the main thrust of the present 
task has been to establish a solution procedure for the scramjet flow. The scramjet flow is 
compressible, turbulent and reacting. The fuel used is hydrogen and the combustion process 
proceeds at a finite rate. As a result, the solution procedure must be capable of addressing such 
flows. The codes chosen for the present work are the TUFF and the SPARK codes. The 
modifications done to the codes (by the principal investigator on SPARK and TUFF and by Dr. 
Gregory Molvik on TUFF) are: 1) added the two-equation k-epsilon turbulence model (with 
low-Re correction terms), 2) added two-equation k-omega turbulence model, 3) added two 
compressibility correction models for high-speed applications, 4) added a 9-species, 20-reaction 
steps finite rate hydrogen-air chemistry model and 5) developed and incorporated a turbulence- 
chemistry interaction model (SPARK only at present). These modifications were carried out 
over the duration of the period mentioned in the beginning. So far, six technical papers have 
been presented (OR to be presented in the near future) in national and international meetings and 
one journal publication has been completed based on the work done under this section of the 
task. Copies of the important publications ( not all of them because the paper presented at two 
international meetings were already published in USA ) are attached at the end of this report. 
Since these publications contain all the relevant technical details of the models, test cases etc. 
those details will not be given in the main body of the report. 



Since the reaction mechanism in the scramjet is highly dependent upon the mixing between the 
fuel and oxidizer streams and the reaction zone is mainly confined (initially) to this mixing 
region, it is logical to validate the solution procedure on reacting mixing layer configurations. 
This is the strategy followed in the present study. First a two-dimensional nonreacting mixing 
layer configuration was computed using all the relevant codes. They produced nearly identical 
results thus proving the basic consistency of the solution procedures. The same geometry was 
used for the mixing and reaction between the fuel and oxidizer streams at high speeds. Since any 
validation procedure must include comparison with experimental data, a search was made to 
secure this data. Useful experimental data in the area of high speed reacting flows is extremely 
rare and available data sets are nowhere near being adequate for a thorough validation effort. 
Available experimental data correspond to two basic configurations: 1) two-dimensional 
(Burrows-Kurkov experiment) and 2) Coaxial jet case where two coflowing jets (inner one is 
hydrogen and the outer one is vitiated air) . Both these cases were used for the validation effort 
as detailed in the attached publications. Representative figures are given in this regard. 

The flow field in the combustor is extremely complex. There are interactions between 
the many different physical aspects of the flow, such as the interaction between the turbulence 
field and the mean flow, interaction between heat release and mean flow and the interactions 
between turbulence and chemistry. Of these, the last one forms the basis for the present work. A 
turbulence-chemistry interactions model was developed in cooperation with the scientists at 
NASA Langley research center (please see attached publications) and this model was 
demonstrated by means of the two-dimensional reacting mixing layer case mentioned in the 
beginning of this section. This model represents the first step in the effort to achieve an accurate 
model for depicting the turbulence-chemistry interactions. The mixing rate is heavily dependent 
upon the turbulence field and if the turbulence field is affected by the chemistry, then it will 
affect the physical aspects of the flow such as ignition (since mixing is affected). Also, from a 
practical point of view, the length of the combustor is a very important parameter in the design of 
the hypersonic vehicle because the propulsion system is highly integrated with the vehicle 
airframe. If the mixing process is affected by the interactions mentioned above, then it could 
lead to impractical dimensions of the combustor. On the other hand, these interactions may have 
a positive effect and help to come up with a more compact geometry for the propulsion system 
than before. The main point here is the fact that the solution procedure must have the capability 
to include these interactions in the computations. The present work provides such a capability 
and the model is being improved. 



The task under the cooperative agreement was designed to last a longer period of time than the 
duration corresponding to this report. As a result, the full tip-to-tail computational demonstration 
on a realistic configuration could not be completed. For example, the combustor flow 
calculations are just beginning and the nozzle flow computations are yet to begin. So this final 
report refers only to the work completed by the principal investigator during the past three years. 
However, the solution strategy is in place and the tools (e.g.. TUFF, grid) are identified and 
readied. Should the need arise, it should be a straight forward process to go from where the 
current work is today to the completion of the full tip-to-tail computations. The principal 
investigator extends the offer to help in initiating such a task should the need arise. 
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1 Comparison of STUFF and UPS forebody solutions (zoom view) 
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FIG. 2 Comparison of STUFF and UPS forebody soujxions 



Forebody Solution 



FIG. 3 Forebody temperature and pressure contours 
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FIG. 5 Grid for hypersonic inlet geometry 






FIG. 6 Grid for hypersonic inlet geometry 
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FIG. 7 2-D reacting mixing layer results 
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ABSTRACT 

The objective of this work is to investigate the ef- 
fect of the interactions between the turbulent scalar 
field and chemistry in compressible, reacting turbu- 
lent flows. A model for these interactions has been 
proposed. In the thermochemistry model, the effects 
of temperature and species concentration fluctuations 
on the species mass production rates are decoupled. 
The effect of temperature fluctuations is modelled 
via a moment model and the effect of concentration 
fluctuations is accounted for using an assumed /?- pdf 
model. A two equation (Ar-u;) model with compress- 
ibility correction is used to calculate the turbulent 
velocity field. Preliminary results for the case of a 
two-dimensional reacting mixing layer are presented. 
Computations are carried out using the Navier-Stokes 
solver SPARK, with a finite rate chemistry model for 
hydrogen-air combustion. 

NOMENCLATURE 

A , b coefficients in Arrhenius rate equation 

Ci A turbulence model constants 

E total internal energy 

f n mass fraction of species n 

g scalar variance (enthalpy, mass fraction) 

H source vector 

H total enthalpy 

h static enthalpy 

k turbulent kinetic energy 

kfjkf, forward and backward reaction rate constants 
L number of reaction steps 

M n Molecular weight of species n 

Mt Turbulent Mach number 

N number of chemical species 

Pr } Pr t laminar and turbulent Prandtl numbers 
p pressure 

R Ci Re x turbulence model coefficients 

Rk,Rui turbulence model coefficients 

Sc,Sc t laminar and turbulent Schmidt numbers 

T temperature 

T a Activation Temperature 

t time 

* Senior Research Scientist, Senior Member AIAA. 


At time step 

U dependent variable vector 

U velocity vector 

w n production rate of species n 

x streamwise coordinate 

Xj j th coordinate 

y transverse coordinate 

a, a* turbulence model coefficients 

ao,c*3 turbulence model coefficients 

6ij Kronecker delta 

€ turbulence energy dissipation rate 

e c compressible dissipation rate 

dissipation rate in ^-equation 
77 compressibility correction coefficient 

7 ratio of specific heats 

specific dissipation rate 
p laminar viscosity 

Pt turbulent viscosity 

v kinematic viscosity 

p density 

Vd^ic, o'* turbulence model constants 
Tj j stress tensor 

$j flux vector in j th direction 

Subscripts 

t turbulent quantity 

INTRODUCTION 

Advanced air-breathing propulsion systems such as 
the supersonic combustion ramjet (scramjet) have 
been studied for a long time as candidates for high- 
speed propulsion. Mixing between fuel and oxidizer, 
heat release, presence of shocks and the chemical re- 
actions are some of the major aspects of such flows. 
Computations and experiments of the flow fields in 
the combustors of such systems reveal a complex in- 
terplay of these physical phenomena. The effect of 
turbulent mixing on chemical reaction is an example 
of such an effect which must be addressed in design- 
ing such a propulsion system. One such effect, the 
turbulence-chemistry interaction, is the subject of our 
present study. Design parameters such as the mode of 
fuel injection, flame stabilization, ignition delay, com- 
bustion efficiency etc. may be influenced by these in- 
teractions. 
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There are many factors that affect the successful 
evaluation of the combustor flow. The flow field in the 
combustor is turbulent and compressible. A further 
complication is introduced by the high heat release 
in the combustor. Today, an ideal procedure for the 
study of such problems (in the absence of exact ana- 
lytical solutions to the governing equations) will be a 
combination of an accurate numerical solver coupled 
with experimental data to (i) validate the numerical 
solver and (ii) to provide valuable comparison data for 
the particular combustor configuration being studied. 
The calculation method must be capable of address- 
ing high-speed, turbulent, reacting, compressible fluid 
flows involving high energy release. There is a glut 
of useful numerical solvers applicable for a wide vari- 
ety of flows including all speed regimes. The capabili- 
ties of the computational fluid dynamics (CFD) codes 
depend upon the sophistication and accuracy of the 
models used to simulate the various physical processes 
involved. There has been a lot of progress made in the 
development of accurate turbulence models in the re- 
cent years. Finite rate chemistry including nonequilib- 
rium thermodynamics is another area where significant 
progress has been made. Computational algorithms 
which are fast and accurate are being improved every- 
day. Computers which can run these codes are also 
being improved. However, every one of these areas 
is still under development and as a result there is no 
single CFD code which includes the accurate models 
and algorithms in all the areas. Limitations such as 
the maximum grid size, temperature range of applica- 
bility of thermodynamics models, number of steps in 
a finite rate reaction model, near- wall applicability of 
the turbulence model etc. still affect the usefulness of 
these CFD codes. On the other hand, experimental 
study of such problems is not possible at present due 
to the lack of adequate facilities. Fortunately, many of 
the major issues can be addressed using current CFD 
tools. 

The focus of the present work is on the turbulence 
model(s) used for reacting, compressible flows. Specif- 
ically, the interaction between the turbulence field and 
the chemical kinetics is of interest. The effect of tur- 
bulence on phenomena such as mixing between two 
streams, boundary layer etc. has been studied exten- 
sively and it is an ongoing process. However, there 
is not enough work done in the area of turbulence- 
chemistry interactions to either emphasize or neglect 
the importance of such interactions. The main reason 
for this is the fact that it is extremely difficult to mea- 
sure these interactions quantitatively. Recently, there 
have been efforts aimed at establishing pdf-based mod- 
els for coupling the turbulence field with the thermo- 
chemical field in numerical simulations. References [1], 
[2], [3], [4] and [5] are just a few of the many reports 


available in this area. There has also been a concerted 
effort to generate validation data bases through direct 
numerical simulations. Even though there is a lot of 
publications in this area, this effort is still in its early 
stages. 

In a computational analysis of flow fields the choice 
of turbulence model(s) is dictated by two conflicting 
requirements. The model should (i) be reasonably ac- 
curate with good physical justification; and, (ii) be 
computationally feasible for complex flow geometries. 
The turbulence model used in the present instance is 
a two-equation (& — u ;) model [6, 7], modified to ac- 
count for compressibility. The effect of compressibil- 
ity in turbulent flows is an important aspect of high 
speed turbulent flows. There has been a significant ef- 
fort aimed at developing models to account for these 
effects [8, 10]. These models have been used to predict 
shear layer growth rates with some success [11, 12]. 
The model used in the present computations is the 
one developed by Zemen [10]. On the thermochemistry 
front, the effects of temperature and species fluctua- 
tions are decoupled. For the temperature-turbulence 
interaction, a moment model is used [1]. The chemical 
species fluctations are accounted for using an assumed 
multivariate-/?-pdf model [1, 2, 3]. The thermochem- 
istry models used here are independent of the model 
for the turbulent velocity field. These models require 
the solution of the evolution equations for all the mean 
flow thermo-chemical variables, turbulent kinetic en- 
ergy (fc), specific dissipation rate (u/), enthalpy vari- 
ance and turbulent scalar energy (sum of all the species 
variances). The CFD code used is SPARK [13], devel- 
oped at NASA Langley Research Center. 

Mixing plays a major role in high speed combustor 
flows. The reaction zone is mainly confined to mixing 
layers that exist between fuel and oxidizer streams. 
Efficient mixing of fuel and oxidizer is very important 
for the design of combustor size. Mixing layers ranging 
from subsonic to supersonic speeds have been studied 
extensively over the years [1, 11, 14, 15, 16, 17, 18, 19, 
20]. However, most of those studies use oversimplified 
models for turbulence-chemistry interactions. In the 
present work, a two-dimensional, compressible, react- 
ing, mixing layer (hydrogen- air) is computed with the 
aforementioned, more realistic turbulence-chemistry 
interaction model. The reaction model (hydrogen-air) 
used is described in Ref. [21] (9 species and 20 reac- 
tion steps) and is given in Table 1. Experimental data 
from compressible, reacting mixing layers is still scarce 
which hinders the validation of the calculation proce- 
dure. As a consequence, the comparison of results is 
confined to those obtained with and without the inter- 
action model. 

The governing and secondary equations used in the 
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computations have all been described in detail in the 
references cited above. Only an abreviated equation 
set will be given in the present paper. The computa- 
tions were performed on the supercomputers of NAS 
and NASA Ames Research Center (C-90). 

GOVERNING EQUATIONS 

The continuity equation, Navier-Stokes equations, 
energy equation and species continuity equations gov- 
ern the instantaneous evolution of the flow vari- 
ables [1, 15, 22], Density- weighted averaging [23] is 
used to derive the mean flow equations from these 
equations. The dependent variables, with the excep- 
tion of density and pressure, are written as 

<t> = $ + 0" (1) 

where the <j>" is the fluctuating component of the vari- 
able under consideration and its Favre-mean <f> is de- 
fined as 

< s > 

In this equation, the overbar indicates conventional 
time-averaging. Density and pressure are split in the 
conventional sense as, 

P = P + P* and p = p + p' (3) 

The averaged continuity and momentum equations are 

VP , Wi n 

at + ~d^~ ~ 0 (4) 


dm dpUjUj dp d/m'/u'f dT~ 

dt dxj dxi dxj + dxj 

where 

, du{ diij 2 dut _ . . 

T “ = ^ + ^-- 3 er t s “ < 6 > 

with repeated indices indicating summation. 

The two turbulence variables for which we carry evo- 
lution equations are the turbulent kinetic energy (&) 
and the specific dissipation rate (w) [6] defined as 

* s PJ w W 


u = e/(C 2 k) 


where 


Rex = 


du" du" 

e = ( 11 ) 

P 

The constant i? c =6. 

Closure of the averaged equations is achieved by 
invoking the Boussinesq approximation which relates 
the turbulent stresses to the mean strain rate. The 
Reynolds stress tensor is written as, 

_ 2 _ 

- pu'i'u'l = fit Sij - - pkSij 

%J dxj dxi 3 dxk t; 

where is the turbulent/eddy viscosity defined as 

fit = ot m p- (13) 


ry* 1 3*7 
a ° + Rk 

1 ^ R h 


ocq — Ci/3, Rk — 8, Ci = — (15) 

The mean continuity equation (4) does not require 
any modeling. The modelled momentum equation is, 


i dpUjUj _ 9p_ 

’ dxj dxi 


, 2 dpk 

3 dxj 


Ct = 


= 9 

100 l + (^f) 4 


There has been a considerable amount of activity in 
the area of modeling the compressibility effects [8, 10]. 
In these models, the compressible dissipation terms are 
expressed as functions of the turbulent kinetic energy 
dissipation rate and the local turbulent Mach number. 
The compressibility effects are represented by a com- 
ponent of the dissipation rate (e c ) given as 

e c = K c e 
K c = rjF(Mt) 

2& 

= -r (it) 

where a is the local speed of sound and F(M t ) is a 
function of the local turbulent Mach number (M t ). As 
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mentioned before, the model used in the present work 
was proposed by Zemen [10] where F(M t ) is given by 

F(M t ) = 1 — exp[—( — Q£*° ) 2 ], M t > M ta 

= 0, M t < M to (18) 

with M to =0.1 and 77=0.75. The modelled turbulent 
kinetic energy equation is (1, 9] 


^ + = *-*.*!+*.> 

+ 1 (19) 


where 


Pk = 


— 77-77 Wi 
pu ‘ u > 


( 20 ) 


The modelled w-equation used in the present anal- 
ysis [6] is given below. This equation does not include 
any compressibility corrections. 


dpw dpcvUj u> n ~ p 

- — + Q x _ = a-j^Pk - Cipuj 2 + <j d ^-D ku 


dt 


where 


and 


k 

d u - fit . dt 

- — [(/i d )— — 

OXj (Tu OXj 


_ dk du> 
k “ ~ dx, dxj 


( 21 ) 


( 22 ) 


Or = 


1 a 0 + 7^ 


2a* 1 + ^ 

ily 


The turbulence constants are 


(23) 


where a $ is a coefficient which, normally, is a constant. 
For <j> = f n (n represents the species) , cr# = Sc t , and 
for the static enthalpy, ( <j> = h) } <j$ = Pr t . 

Using the above definition, and omitting the body 
force contribution, the time-averaged and modelled en- 
ergy equation [1] is 


dpE dpEU j d , 

~dT + = {Tij ~ p 6i > ~ p u >'D u ' 

, <L\tA , h \ ** — tor', 

dxj [( Pr + Pr t ) dxj + ^ + a J d Xj * 26) 


where a k is a coefficient that appears in the turbulent 
kinetic energy equation. The modeled species continu- 
ity equation is 


dpfn dpfn Uj 

dt dxj 




In the above, the mean species production rate due 
to chemical reaction (u> n ) needs to be modelled. This 
term is a function of both the temperature and species 
concentrations. In the absence of a model for the 
interaction between the turbulent temperature and 
concentration fields, this term is evaluated using the 
mean values of local temperature and species concen- 
trations. In a finite-rate system involving L reaction 
steps and N species, the instantaneous production rate 
of a species n can be represented (law of mass action) 
in the following general form: 


l On 


L 
;=i 


(28) 


N 


(‘w”' -w"' ntTf)-"')' 


N 


f. 


i=i 


M, 


M, 


a 0 = 0.1, 
Pr = 0.72, 

and 


Ru, = 2.2, 

Pr t = 1.0, 


<?k = 1-0, 

Sc = 0.22 


<r u = 1.67, 
Sc, = 1.0 


where 


N 


N 


m, 


= E ni = E v > 




5 = 1 


<Jd = 0 . 0 , D ku < 0 

= 0.3, D kw > 0 


The mass-averaged total energy can be written in 
terms of the total enthalpy as 


E = H - l (24) 

P 

The correlations between the fluctuating velocity and 
the scalar fluctuations are modelled using a gradient- 
diffusion hypothesis. A typical model is of the form 


— pu , i ( <j) ,f 


El (EE) 

<T<t> dxi 


(25) 


In the above equations, and i/'j are the number of 
molecules of the scalar s involved in the /-th reaction 
step in the forward and backward directions, respec- 
tively. The forward and backward rate-constants of 
the reaction / are given by k/i and kbi respectively. 
The reaction rates are usually strong functions of the 
temperature: 


k Jt = A,T b < exp{ --f] (29) 

where A/, b\ and T ai are numerical constants specific 
to the given reaction step /. 

Turbulence- chemistry interaction model 
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The interactions between the temperature and con- 
centration fields (turbulent) may be of importance in 
many applications. The following section describes a 
model for such interactions. Practical combustor flows 
involve multiple scalar mixing and reactions. In order 
to calculate the mean species production rate using 
the pdf approach, a model for the joint pdf of temper- 
ature and various species is required. Such a joint pdf 
is difficult to specify and as a preliminary step we effect 
the simplification that temperature and species fluctu- 
ations are uncorrelated. This permits us to deal with 
the temperature and species fluctuations separately. 

Equation (29) can be written as 


k fl = A,(f + T”) b ' exp[- 


c T + T") 


(30) 


Assuming that XL < 1, the term (T+T") can be ex- 
panded in a series and the resultant modified reaction 
rate term is written as, 


k h = (1 + m) A,T b ‘ exp[— S-] (31) 

where 


m 


[(6,_1) ( b j + ^) + 


17^, 2 

2 { f f 2 


(32) 


Terms of order higher than two in -C- are neglected 
from the series expansion in the present analysis since 
accurate models for these higher order terms do not ex- 
ist. The factor m represents the effect of temperature 
fluctuations on the mean reaction rate. The tempera- 
ture variance is calculated from 


’J'n r £n 


h"h" 



(33) 


Decoupling the effects of temperature and species 
concentrations, the mean species production rate can 
be written as, 

L 

W n = M n — V* n i) X (36) 

(=1 

_ N __ N 

{ vdl - vdl M r“‘)h'} 

5 = 1 5=1 

where 

5 = 1 5 = 1 

In the above, angular brackets represent conventional 
time averaging. The modelled expression for Ifi using 
the multivariate /?-pdf is [2], 

N v ,i mi 

i i‘ - n ii(^* + "w - r >/ ri( 5 + m ' ~p) ( 3? ) 

5 = 1 r=l p = l 

Similarly, the expression for hi is, 


N v"x n x 

ibt - n m + v "i - r )/ n<* +*»-*> 

5 = 1 r=l j>=i 

where 

B = f3 1 + ft + • • ■ + &n 


(38) 

(39) 


The production of turbulent scalar energy due to 
chemical reaction can be decomposed as 


iV S N 

E Wn K = ^Wnfn ( 4 °) 

n = l n = l n = l 

Only the first term on the right hand side of the above 
equation needs further modeling. This term can be 
written as, using the assumptions above, 


in the present work. Other pdf-based models for the 
temperature effects exist [5] which will be explored in 
the future. 


For reacting flows involving multiple chemical 
species, Girimaji suggests [2, 3] the use of a multivari- 
ate /?-pdf model to account for the effects of the scalar 
fluctuations on the species production rates. This 
model is briefly outlined below. The parameters of 
the multivariate /?-pdf for the iV-scalar mixing process 
(A > * * * , 0n) are functions of the mean mass fractions 
f n and turbulent scalar energy Q : 


where 


2* 

II 

>1 


(34) 

N 

N 


E 7 » 2 . 

n = l 

Q=T, fn fn ' 

n = l 

(35) 


N AT L 

E = E E^"' - V nl) x ( 41 ) 

n = l n = 1 (=1 

__ JV iV 

{k„p m '(Y[Mr->)j :u ~ w»'(f[ 

5=1 5=1 

where 

J /U = (fn n /* , ")> = (fn n f*"') 

5=1 5=1 

The terms Jfi n and are also moments of the scalar 
joint-pdf. The modelled expression for these moments 
are [2], 


J/U = Iji 


Jbi„ — hi 


Pn + v' nl 
B + m ; 
Pn+jC 

B + Til 


(42) 


5 



Substitution of equation (42) into equation (41) leads 
to the model for the source/sink of turbulent scalar 
energy. The main advantage of this choice of the 
assumed-pdf model is that the chemistry related mod- 
els are obtained analytically and no numerical inter- 
gration in the species space is necessary. 


The turbulence-thermochemistry interaction models 
require the enthalpy variance and the turbulent scalar 
energy distributions. The modeled equations for en- 
thalpy variance (h"h") and turbulent scalar energy(Q) 
are of a form similar to that of the turbulent kinetic 
energy: 


dpg 

dt 


+ 


+ 


dpgUj 

dxj 


- 2 pu''G 


dG 

dxj 


)c:K5 + £>& + # 


dxj 


dxj 


2pc g 


(43) 


For g = G = h, #=0, cr — Pr and <r g = Pr t . 

For 9 = Ell fUn, G = /n, 2El t <7 = 
Sc and <x g = Sc t . The dissipation term in the above 
equation is assumed to be 


*g ~ Cg^9 = C g C 2 ug (44) 

The model for ™n/n » is given in equations (40)- 

(42). The constants C g are assigned a value of 0.5. 


Solution of the modeled governing equations 


The equations are discretized and integrated in 
space and time to obtain steady state solutions using 
an elliptic solver SPARK [13]. The governing equa- 
tions are written in vector form as follows. 


dV d$j 
dt dxj 


= H 


(45) 


where U is the vector of dependent variables, are 
flux vectors containing convective and diffusive terms 
(repeated indices indicate summation), and H is the 
source vector containing production/dissipation terms. 
The temporally discrete form of equation (45) is 


d 

U n+1 = U" - A t [-£-*- - H” +1 ] (46) 

OXj 

where n is the old time level and n + 1 is the new time 
level. 


The source terms in the k- and ^-equations are de- 
coupled by suitable manipulation of the lj term in the 
present analysis in order to alleviate the computational 
stiffness introduced by these terms. For example, in 
the ^-equation, the dissipation term is written as, 


C 2 pwk — C 2 pk u* 


(47) 


The term c j* is taken from the most recent calcula- 
tion step. The source terms of the w-equation are 
also manipulated in a similar manner. These nonlin- 
ear turbulence source terms are treated in a pointwise- 
implicit manner while solving the turbulence equations 
by rewriting equation (46) as, 

(!- (U n+1 -U n ) = - At - H n ] (48) 

The discretized equations are solved by means of a 
fourth-order compact scheme. 

RESULTS AND DISCUSSION 

Choosing a flow configuration to demonstrate the ef- 
fect of turbulence-chemistry interactions is a difficult 
task due to the lack of guidelines based on prior data 
(experimental or otherwise). Since the present work 
is part of a larger task of establishing a solution pro- 
cedure for high speed propulsion systems (scramjets), 
the initial choice fell on a two-dimensional high speed 
reacting mixing layer. Non-reacting and reacting high 
speed mixing layers have already been computed us- 
ing the computational procedure described above [14]. 
Available experimental data were used to validate the 
prediction procedure. The /?- pdf model to account for 
the effect of species concentration fluctuations on the 
mean production rate of species was also introduced in 
conjunction with the k-e turbulence model before [1]. 
The main aim of the present work is to introduce the 
results obtained by coupling the /?-pdf model with an 
improved version of the k — u model [6]. 

A schematic of the flow problem is given in Figure 1. 
The two streams are, air ( U =1606 m/sec, T =1600 K 
with /*, =0.0, / 0a = 0.267 and /*, =0.733) and hy- 
drogen (U =1250 m/sec, T =254 K with /*r, = 1.0, 
/o 3 = 0.0 and /at 3 = 0.0). The two streams are super- 
sonic with the air stream Mach number of 2.07 and 
hydrogen stream Mach number of 1 .03. The inlet mean 
velocity is assumed to have a hyperbolic tangent pro- 
file, thus imitating the flow that exists downstream of 
the splitter plate trailing edge. A constant turbulence 
intensity level is used in the free stream for arriving 
at the initial distribution of turbulent kinetic energy 
and the specific dissipation rate. The pressures are 
matched between the two streams (P =1 atm.). A 13- 
step, 8-species H 2 — Air reaction model (Table 1) has 
been used for the finite-rate chemistry system consid- 
ered here. A 101 X 81 grid (101 points in flow direc- 
tion, 81 points in the transverse direction) was used 
for the calculations. The length of the flow domain 
is 0.25 m and its width is 0.05 m. In all the figures 
shown in this report, y refers to the lateral distance 
measured from the outer edge of the lower stream. The 
title g refers to calculations that include the interac- 
tion model. 
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As a first step, the difference made by the temper- 
ature fluctuations model ( 28) is explored. Figure 2 
shows the factor (1-f-m) as a function of y for repre- 
sentative reaction steps at the exit plane of the flow 
domain. The need for including the temperature fluc- 
tuations is evident in the figure. The forward reac- 
tion rate changes by as much as 400 percent (reaction 
step 1) due to the effect of this model. The species 
production rate is a strong function of the reaction 
rate ( 32) and hence will be affected by the inclusion 
of this model. 

The effect of the turbulence-chemistry interactions 
model on the species production rate is shown in the 
next few figures. Figure 3 shows the production rates 
of the major species involved (Ho, 0 2 , HoO , OH) 
as functions of y at the exit plane. Computational 
results obtained with (k-w-g) and without (k-w) the 
effect of these interactions are shown in the figure. It 
is seen that the net change in the production rate in 
the shear layer/reaction zone is reduced by the effect 
of these interactions for the species H 2 , O 2 and HoO . 
For the species OH, the effect of these interactions 
seems to be to reverse the rate of production from the 
case where the interactions are not included. While 
it is not conclusive whether the turbulence-chemistry 
interactions increase or decrease the production rates 
of individual species, it is important to note that these 
interactions do have a significant effect on progress of 
the chemical reactions. 

Figures 4 show the distribution of the streamwise 
velocity and static temperature as a function of y at 
the exit plane. The effect of the interactions seems to 
be felt more at the edges of the mixing layer/reaction 
zone than anywhere else. In the Going back to fig- 
ure 3, the production rate seems to be sensitive to the 
temperature gradients rather than the temperature it- 
self. Turbulence affects the distribution of tempera- 
ture both directly (thermodynamic energy equation) 
and indirectly (via the shear layer effect). This effect 
is transmitted through the interactions model and is 
seen in the chemical production rate profiles. Figure 5 
shows the effect of the interactions on the species dis- 
tributions. The effect seems to be more pronounced in 
the case of the products (H 2 0, OH) than the primary 
reactants (H 2 , 0 2 ). The changes are larger in regions 
of higher gradients (on the fuel stream side for H 2 and 
air stream side for O 2 , for example). 

One of the unknown yet crucial factors associated 
with the present work is the lack of prior knowledge 
as to the effect of turbulence on chemical reactions 
and the reverse effect of the chemical reactons (heat 
release, concentration change) on the turbulent field. 
Which of these effects is dominant is a debatable is- 
sue. Based on the results shown in this report as well 


as other similar computations (not shown) carried out 
by the author, it is not possible to come up with a def- 
inite answer to this question. It is also not established 
whether the type of flow configuration (mixing layer, 
jet, boundary layer) has any significance in this dis- 
cussion. So far, all computations using the proposed 
model have been done for reacting mixing layer (2- 
D and coaxial jets) involving the mixing and reaction 
between fuel and oxidizer streams. It is not certain 
whether a premixed reacting flow configuration will 
show a amore pronounced effect of these interactions 
than the non-premixed cases studied so far. Also the 
effect of the interactions near a solid wall is an area 
worth exploring. 

An interesting aspect of the analysis is that the par- 
ticular turbulence model chosen (k - e as opposed to 
k — w, for example) seems to have an effect on the 
overall predictions [1], Since the ^-equations, which 
provide the variances needed to construct the interac- 
tions model, are strongly coupled (via their produc- 
tion terms) with the turbulence model equations it is 
also essential to differentiate between the effects of the 
turbulence model and the turbulence-chemistry inter- 
action model. This brings up the importance of the 
accuracy of the turbulence model which was one of 
the main reasons for carrying out this task using the 
k — <jj turbulence model. It is crucial to keep the short- 
comings of the model in perspective while attempting 
to validate the results. The assumptions made in ar- 
riving at the model, such as the decoupling of temper- 
ature and concentration effects, are as important as 
the model itself in some cases. 

On the positive side, this effort represents a signifi- 
cant push in the right direction in this very important 
area of turbulent reacting flows. The model is simple 
to use and is easily adaptable to other types of tur- 
bulence closures such as the Reynolds stress models. 
More improvements can certainly be done in terms of 
using more realistic models for the temperature and 
concentration fluctuations such as a pdf-based model 
for temperature effects (as opposed to the moment 
model employed here), a joint-pdf model to represent 
the effects of both the temperature and the concen- 
tration fields etc. However, in order for the model to 
be useful for practical applications, it must be kept in 
a form which promotes ease of use as well as compu- 
tational economy. This is the main reason for resort- 
ing to a two-equation level turbulence model in the 
present work. A concerted validation effort is crucial 
to ensure the success of the modelling efforts. One 
debilitating factor in this regard is the nonexistence 
of useful experimental data. It is extremely difficult, 
if not almost impossible, to obtain such experimen- 
tal data using present day equipment. An alternative 
seems to be direct numerical simulations (DNS) which 
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has shown promising signs in relatively simpler flow 
configurations. Inspite of the major advances made 
in the area of DNS in recent years, it still is a devel- 
oping field and is far from providing useful data for 
validation purposes in cases such as the present work. 

CONCLUSIONS 

A turbulence-chemistry interaction model has been 
proposed in conjunction with a recent version of the 
two-equation k — w model of turbulence for use in 
chemically reacting flows. Preliminary computations 
carried out with the model for the case of a two- 
dimesional, high speed, reacting mixing layer indicate 
that these interactions have a significant effect on the 
flow predictions. The species production rates of in- 
dividual species seem to be affected by these interac- 
tions. It is not possible, yet, to quantify the effects 
of these interactions based on available data. More 
detailed analysis of the various aspects of the model 
developmental process need to be done. 
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Table 1. H 2 — Air Reaction System 


No. 

Reaction 

1 

H + 0 2 — O + OH 

2 

OH + H 2 ~ H 2 0 + H 

3 

o + h 2 — oh + h 

4 

OH + OH — H 2 0 + 0 

5 

H + OH + M“ H 2 0 + M 

6 

h + h + m*=h 2 + m 

7 

o + o + m^o 2 + m 

8 

H + O + M^OH + M 

9 

H + 0 2 + M ** H0 2 + M 

10 

OH + H0 2 ** h 2 0 + 0 2 

11 

H + H0 2 *±H 2 + 0 2 

12 

h + ho 2 -oh + oh 

13 

0 + ho 2 * oh + o 2 

14 

ho 2 + ho 2 h 2 o 2 + 0 2 

15 

H + H 2 0 2 -H 2 + H0 2 

16 

OH + H 2 0 2 ** H 2 0 4- H0 2 

17 

H + H 2 0 2 - H 2 0 + OH 

18 

O + H 2 0 2 m H0 2 + OH 

19 

OH + OH + M ** H 2 0 2 + M 

20 

OH + OH m H 2 + 0 2 


Species : H 2f 0 2i H 2 0 } OH , H, O, H0 2) H 2 0 2 and 
A^inert) 

M is a third body (all species included) 



Fig. 1 Mixing Layer - schematic 
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ABSTRACT 

Computations are done on flow configurations that 
resemble the reaction zone in the scramjet combus- 
tor flows. Compressible, reacting, turbulent flow so- 
lutions are obtained. A two equation (fc-e) model 
with compressibility correction is used to calculate the 
flow field. A finite rate (8-species, 13-reaction steps) 
chemistry model for hydrogen-air combustion has been 
used. Computations are carried out using the Navier- 
Stokes solver TUFF. Predictions are compared with 
available experimental data and also those obtained 
by using the code UPS. 

NOMENCLATURE 

A, 6 coefficients in Arrhenius rate equation 

Ci,C 2 ,CV turbulence model constants 
E total internal energy 

f n mass fraction of species n 

H total enthalpy 

h static enthalpy 

k turbulent kinetic energy 

k/.ki, forward and backward reaction rate constants 
L number of reaction steps 

M n Molecular weight of species n 

M t Turbulent Mach number 

N number of chemical species 

Pr,Pr t laminar and turbulent Prandtl numbers 
p pressure 

Sc,Sc t laminar and turbulent Schmidt numbers 
T temperature 

T a Activation Temperature 

t time 

U velocity vector 

w n production rate of species n 

x streamwise coordinate 

* Senior Research Scientist, Senior Member AIAA. 

tPostDoctorla Fellow, Member AIAA. 

♦Senior Research Scientist, Member AIAA. 


xj j th coordinate 

y transverse coordinate 

Sij Kronecker delta 

e turbulence energy dissipation rate 

e c compressible dissipation rate 

rj compressibility correction coefficient 

7 ratio of specific heats 

uj specific dissipation rate 

/i laminar viscosity 

Pt turbulent viscosity 

v kinematic viscosity 

p density 

<7 e turbulence model constants 
T{j stress tensor 

flux vector in j th direction 

Subscripts 

t turbulent quantity 

INTRODUCTION 

Hypersonic travel requires propulsion systems which 
are different from the conventional ones used in most 
of the modern aircraft. The supersonic combustion 
ramjet (scramjet) is a system considered to be suit- 
able for high speed applications. There has been a 
tremendous amount of activity in the area of scramjet 
research in recent years ( [1] - [11]). Some of the re- 
lated topics include inlet configuration, mixing layers, 
mixing enhancement, combustor configuration, finite 
rate chemistry models and chemical kinetics. The fuel 
used in the scramjet varies depending upon the ap- 
plication. For example, hypersonic waveriders using 
hydrocarbon fuels have been designed [12] for applica- 
tions in the moderate hypersonic speed regimes. For 
Mach numbers of the order of 15 and above, hydro- 
gen is generally considered to be the fuel of choice. 
In the present work, hydrogen is the fuel used in the 
computations. 

The present work represents a computational ef- 
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fort in establishing a solution procedure for hypersonic 
propulsion applications. The entire task of establish- 
ing the solutions procedure must then be divided into 
smaller tasks dealing with subsets such as turbulence 
modelling, chemical kinetics, geometry etc. One such 
task is the topic for the present study. Here, the rele- 
vant flow features of the combustor, namely the mixing 
and chemical reaction between the fuel and oxidizer 
streams, is addressed. The flow field in a scramjet is 
complex. It is turbulent and compressible involving 
high heat release. The solution procedure should ad- 
dress all aspects of the flow field adequately. It should 
be capable of accurately modelling the turbulent field, 
taking into account the effects of compressibility, and 
addressing the changes associated with heat release. 
Also, the interactions between the distinct physical 
aspects of the flow such as the effect of heat release 
on turbulence, the interaction between turbulence and 
chemistry etc. must be properly addressed. Signifi- 
cant progress has been made in addressing these areas 
via accurate and realistic modelling in recent years [6]. 

Remarkable advances have been made in the area 
of turbulence modelling, accounting for a variety of 
factors that affect the flow field. Compressibility cor- 
rection models to account for the effects of compress- 
ibility, near-wall turbulence models to deal with the 
transition from fully turbulent to zero turbulence, vis- 
cous dominated flow field near no-slip boundaries and 
modifications to models to account for flow curvature 
are some examples. A wide variety of turbulence mod- 
els, including algebraic (zero-equation), one-equation, 
two-equation, Reynolds stress and large eddy simula- 
tion models, are available (for example references [13] - 
[15]) depending upon the sophistication and accuracy 
desired and the limits imposed by numerical solution 
procedures. 

Thermodynamic and chemical kinetic models [16] 
applicable to the scramjet flows have been undergoing 
continuous improvements in recent years. Accurate 
modelling of thermodynamic variables as functions of 
temperature which are valid over a wide range of tem- 
peratures is an example. In flows such as the one 
associated with the scramjet, the time scales associ- 
ated with fluid dynamics and chemical reaction (not to 
mention the turbulence scales) require that the com- 
bustion process be modelled via a finite rate chem- 
istry mechanism. Such a mechanism should account 
not only for the major species (reactants and prod- 
ucts) involved in the chemical reactions but also the 
intermediate transient ones which play a vital role in 
the reaction progress process. Accurate models for the 
chemical reactions in the scramjet combustor, thus, is 
a crucial aspect of the solution procedure. 

The design of the combustor is strongly dependent 


upon factors such as the mixing between fuel and ox- 
idizer streams, presence of shocks in the flow field, 
boundary layer effects, flow separation, extent of chem- 
ical reaction within the combustor and so on. The nu- 
merical solution procedure should have the capability 
of addressing all of these factors while maintaining the 
required accuracy and robustness. There is a glut of 
useful numerical solvers applicable for a wide variety 
of flows including all speed regimes. Computational 
algorithms which are fast and accurate are being im- 
proved everyday. 

Even though there are a wide variety of sophisti- 
cated and physically accurate thermodynamic, chem- 
ical kinetic and turbulence models available it is not 
always possible to use the most accurate and elabo- 
rate versions in a numerical simulation due to the lim- 
itations imposed by computer memory requirements, 
computational economy, ease of use and adaptability 
to practical problems. Solutions often are required, 
especially in ihe engineering industry which is the end 
user for such solvers, in a short time using comput- 
ers that may not be the fastest available. As a result, 
compromises must be struck between physical accu- 
racy and computational feasibility and it is this aspect 
which differentiates between various solvers that exist 
today. 

In the present study, an attempt is made to es- 
tablish a solution procedure for scramjet combustor 
flow predictions from the perspective of the discus- 
sion above. The models chosen to represent the tur- 
bulent and chemistry fields reflect the compromise be- 
tween physical accuracy and computational economy 
mentioned above. The code chosen for the computa- 
tions is the TUFF [17] code and the solutions are com- 
pared with those obtained with the UPS [18, 20] code. 
The turbulence model chosen is the two-equation k — e 
turbulence model with low Reynolds number modifi- 
cations [13]. However, the Baldwin- Lomax algebraic 
model is also available as an option. The compress- 
ibility effects are included via the compressibility cor- 
rection model proposed by Zeman [19]. The fuel used 
is hydrogen although the numerical solver can easily 
be modified for hydrocarbon fuels. A 9-species, 20- 
reaction steps chemical kinetics model for hydrogen- 
air combustion [16] is available. For the computations 
presented in this report, an abbreviated version (8- 
species, 13-steps) of this model has been used. 

Mixing plays a major role in high speed combus- 
tor flows. The reaction zone is mainly confined to 
mixing layers that exist between fuel and oxidizer 
streams. Two flow configurations are chosen for the 
study. The first is the well known Burrows-Kurkov 
experiment [21] in which hydrogen and vitiated air 
streams (two-dimensional) mix and react. The second 
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case is that of an axisymmetric configuration [22, 23] 
where two coaxial jets (fuel and oxidizer) mix and re- 
act. Experimental data from compressible, reacting 
mixing layers is still scarce which hinders the valida- 
tion of the calculation procedure. The available data 
from the above two experiments are used to compare 
with the predictions. The governing and secondary 
equations used in the computations have all been de- 
scribed in detail in the references cited above. Only an 
abreviated equation set will be given in the present pa- 
per. The computations were performed on the super- 
computers of NAS and NASA Ames Research Center 
(C-90). 


GOVERNING EQUATIONS 

The equations used for computations are described 
in detail in references [6, 7, 10] and [24]. Only the 
forms of the modelled equations used in the present 
study are given here. Density-weighted averaging is 
used to derive the mean flow equations from the in- 
stantaneous coservation equations. The dependent 
variables, with the exception of density and pressure, 
are written as 

0 = 0-f<£" (1) 

where the 0" is the fluctuating component of the vari- 
able under consideration and its Favre-mean 0 is de- 
fined as 

? = 4r (2) 

P 

In this equation, the overbar indicates conventional 
time-averaging. Density and pressure are split in the 
conventional sense as, 


and 


€ 


du" du " 


( 8 ) 


Boussinesq approximation is used to obtain closure 
of the averaged equations. Here the Reynolds stress 
tensor is written as, 


Sij 


piti itj — Sij ^ pkSij 

dUj dUj_ 2 dUl 

dxj dxi 3 dxk 


(9) 


where p t is the turbulent/eddy viscosity defined as 


fit = C>— (10) 

with Cn= 0.09. 

The modelled momentum equation, then, is, 

dpUj dpUjUj _ _ dp_ 2 dpk 

dt dxj ~ dxi + 3 &ij 

“ Q x . [ A) Sij] (11) 


The effects of compressibility are included via the 
model proposed by Zeman [19]. Here, the compress- 
ible dissipation terms are expressed as functions of the 
turbulent kinetic energy dissipation rate and the local 
turbulent Mach number. The compressibility effects 
are represented by a component of the dissipation rate 
(e c ) given as 


p = p + p f and p = p + p f (3) 

The averaged continuity and momentum equations are 


d? dpUi 

dt dx{ 


(4) 


dpUi dpUjUj _ dp_ 

dt + dxj ” dxi 

where 






dpu'/uj 

dxj 


% »> 


- X 

3 dx k 6ii 


( 6 ) 


with repeated indices indicating summation. 

In the two-equation turbulence model, the two tur- 
bulence variables are the turbulent kinetic energy (fc) 
and the dissipation rate (e) [13] defined as 


k = 

2~p 


(7) 


(c = K c e 


K c = tjF(M,) 



( 12 ) 


where a is the local speed of sound and F(M C ) is a 
function of the local turbulent Mach number (Mt). 
F(M t ) is given by 


F(M t ) = l-exp[-( \^ n M t > M to 
= 0, M t <M t0 

with A/ to =0.1 and 77=0.75. The modelled turbulent 
kinetic energy equation is [6, 13] 


dpk dpkilj 
dt + dxj 


-pc(l + K e ) 

:!» + *>£ 


(13) 
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where 


where 


p -77-77 dU i 

p> = ~ “/ 8^ 


( 14 ) 


The modelled e-equation used (no compressibility 
corrections) in the present analysis [13] is given below. 


dpe dpeUj 
dt + dxj 


+ 


(Cx P k - C 2 pe) - 


where P k is the production term in the turbulent 
kinetic energy equation. The model constants used 
in the analysis are Ci=1.44, 072=1.92, < 7 *= 1.0, 

^€=1-3, Pr= 0.72, Pr,=1.0, Sc=0.22and Sc, =1.0. 

The mass-averaged total energy can be written in 
terms of the total enthalpy as 


N y 

m ' = n, = 

1 3 = 1 

where, v\ x and t/” ( are the number of molecules of the 
scalar s involved in the /- th reaction step in the forward 
and backward directions, respectively. The forward 
and backward rate-constants of the reaction l are given 
by kji and k^i respectively. 

h‘ = A,T b < exp[ -Tf] (21) 

where A/, b/ and T a , are numerical constants specific 
to the given reaction step /. kn is determined from the 
equilibrium constant for the /-th reaction step and k !t . 

Solution of the modeled equations 


E = H- V - (16) 

p 

The correlations between the fluctuating velocity and 
the scalar fluctuations are modelled using a gradient- 
diffusion hypothesis. A typical model is of the form 


- P u'l4>" 


ii ( QL) 

<7* dxi 


(17) 


where is a coefficient which, normally, is a constant. 
For <j> = /„ (n represents the species) , = Sc t , and 

for the static enthalpy, (<f> = h), a $ = Pr t . Using the 
above definition, and omitting the body force contri- 
bution, the time-averaged and modelled energy equa- 
tion [6] is 


dpE , dpEUi d . 

~df + (T,J “ p ~ rf’vjM 

, 9 u H 1 \ dh ji t dk 

d Xj [ (pr + Pr t } dx } +iP+ a k } d^ ll8) 


where a k comes from the turbulent kinetic energy 
equation. The modeled species continuity equation is 


dpfn dpf n Uj 

dt dxj 





The modelled form of the mean species production 
rate due to chemical reaction ( w n ) is given, for a finite- 
rate system involving L reaction steps and N species, 
in the following general form: 


= M n X>", - *',) X 

1=1 


( 20 ) 




The equations are discretized and integrated in 
space and time to obtain steady state solutions using 
the finite- volume based numerical solver TUFF [17]. 
The TUFF code contains many desirable features for 
the computation of three-dimensional, hypersonic flow 
fields. It has non-equilibrium, equilibrium and perfect 
gas capabilities along with an incompressible option. 
It employs a finite-volume philosophy to ensure that 
the schemes are fully conservative. The upwind in- 
viscid fluxes are obtained by employing a new tem- 
poral Riemann solver that fully accounts for the gas 
model used. This property allows the flow field dis- 
continuities such as shocks and contact surfaces to be 
captured by the numerical scheme without smearing. 
Total Variation Diminishing (TVD) techniques are in- 
cluded to allow extension of the schemes to higher or- 
ders of accuracy without introducing spurious oscilla- 
tions. The schemes employ a strong coupling between 
the fluid dynamic and species conservation equations 
and are made fully implicit to eliminate the step-size 
restriction of explicit schemes. This is necessary since 
step-sizes in a viscous, chemically reacting calculation 
can be excessively small for an explicit scheme, and 
the resulting computer times prohibitively large. A 
fully conservative zonal scheme has been implemented 
to allow solutions of very complex problems. The 
schemes are made implicit by fully linearizing all of 
the fluxes and source terms and by employing a mod- 
ified Newton iteration to eliminate any linearization 
and approximate factorization errors that might oc- 
cur. Approximate factorization is then employed to 
avoid solving many enormous banded matrices. As 
mentioned before, the options for turbulence models 
include both zero and two equation models (both k - c 
and k - w). For more details about the solution pro- 
cedure the reader is directed to the reference cited 
above [17]. 
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RESULTS AND DISCUSSION 

Two reacting flow configurations have been chosen 
for the present study. As mentioned before, the Navier 
Stokes solver TUFF has been used for the computa- 
tions. The first one is the case of coaxial jets [22, 23] 
where a hydrogen jet flows (inner jet) coaxially with 
an outer vitiated air (mass fractions: oxygen=0.246, 
water=0.209 and nitrogen=0.545) jet. A schematic of 
the flow problem is given in Figure 1. The two streams 
are, air (t/= 1380 m/sec, T=1180 K with p=107000 
N/m2) and hydrogen (U =1774 m/sec, T — 545 K with 
p — 112000 N/m2)., The air stream is supersonic with a 
Mach number of 1.97 and the hydrogen stream Mach 
number is LOO. The inlet mean velocity is assumed 
to have a step profile with the two jets having uni- 
form speeds at the specified values (no experimental 
data available). The velocity in the lip region of the 
inner jet tube wall (finite wall thickness) is assumed 
to be zero. The inlet temperature profile is derived 
based on the experimental data given for a location 
just downstream of it (shown later). The inlet species 
mass fraction distributions are also chosen based on 
the experimental data provided at the same down- 
stream location. A constant turbulence intensity level 
is used in the free stream for arriving at the initial 
distribution of turbulent kinetic energy and the dis- 
sipation rate. A 13-step, 8-species Hi — Air reaction 
model (Table 1) has been used for the finite-rate chem- 
istry system considered here. A 81 X 91 grid (81 points 
in flow direction, 91 points in the radial direction) was 
used for the calculations. The inner jet/tube diameter 
(D=0. 00236 m) is used as a reference length. The to- 
tal length of the flow domain is equal to 43.1 D . The 
outer boundary (radial) of the flow domain is taken to 
be at y—ll D . A more detailed description of the flow 
parameters is given in Table 2. The region outside the 
limits of the air jet is assumed to be still air at a tem- 
perature of 273 K. The two-equation (k — c) turbulence 
model is used along with the finite rate H 2 -Air chem- 
istry model mentioned above. In all the figures shown 
in this report, y refers to the radial distance measured 
from the axis of the coaxial system of jets. 

Figures 2-3 show the results of the computations. 
Figure 2 shows the computed and experimental distri- 
butions of species mole fractions. The figure is de- 
signed in a two-column format. The left side col- 
umn represents the inlet (first x-location) data and 
the right side column is the data at the exit plane 
(x/D= 43.1 D). As seen in these figures, the inlet 
data agreement between the computations and exper- 
iment is not perfect, especially around the jet edges, 
and this might affect the computed distributions at 
downstream locations. The comparison between pre- 
dictions and experiment at the downstream location 


(x/D=43.1 D) is good given the above mismatch be- 
tween the two data at the inlet. The development 
of the reaction zone after ignition is not predicted 
well. The experimental data indicates that the reac- 
tion zone (depicted by the water mole fraction distri- 
bution) spreads more quickly than the predictions in- 
dicate. The predictions show the reaction zone to be 
off-center whereas the experimental data shows the re- 
cation zone to be closer to the axis of symmetry. How- 
ever, there is very good qualitative agreement berween 
the data with the peak values of the reaction prod- 
ucts predicted very well. The flow domain was seen 
to have a wave-like structure as shown by the pre- 
dicted profiles. The worst agreement seems to be for 
the case of oxygen. However, when the initial pro- 
files of oxygen are compared one finds that there too 
is the worst agreement between computations and ex- 
periment which may be reason for the problem down- 
stream. Figure 3 shows the comparison of static tem- 
perature data. The agreement between predictions 
and experiment is good qualitatively displaying similar 
trends. The uncertainty associated with the accuracy 
of the experimental data is unknown. There are con- 
siderable differences between the data presented by the 
two references [22, 23], especially in the temperature 
profiles. Overall, there is good qualitative agreement 
between the predictions and experiment. 

The second test case considered is the Burrows- 
Kurkov experiment [21]. The flow configuration is two- 
dimensional. A schematic diagram of the configuration 
is given in figure 4. No-slip walls bound both the upper 
and lower regions (y=0 and y=y maar ). The lower wall 
is inclined to form an expansion surface. Hydrogen is 
injected along this surface into a vitiated air stream. 
The two streams mix and react downstream of the 
injection location (inlet). The hydrogen stream is in- 
jected at a velocity of 1216 m/sec and a temperature 
of 254 K. The airstream comes in at a speed of 1764 
m/sec and a temperature of 1270 K. Full details about 
the flow parameters and geometry are given in Table 3. 
In this case, the reference length used in the hydrogen 
jet width at inlet, h (/i=0.004 m). The models for tur- 
bulence and chemistry are identical to the ones used 
for the coaxial jet case. The grid size is 81 X 121 
(81 grid points in the axial (x) direction and 121 grid 
points in the transverse direction). The total length 
of the solution domain is 0.356 m (x//i=89). Avail- 
able inlet data have been used for the first-plane pro- 
files which improved the predictions remarkably over 
the solutions obtained with uniform profiles. The so- 
lutions are compared with the available experimental 
data at this location (exit) in figures 5-7. The solu- 
tions carried out with the space marching PNS code 
UPS [20] using the Baldwin-Lomax turbulence model 
are also given for comparison. 
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Figure 5 shows the compariosn between the pre- 
dicted distributions of the species mole fractions and 
the corresponding experimental data. As seen in 
these figures, there is excellent agreement between the 
TUFF predictions and experiment. The predictions by 
the UPS code do not agree very well but still there is 
very good qualitative agreement with the experimental 
data. Figure 6 compares the predicted profiles of exit 
plane total temperature and Mach number with the 
experimental data. There is good qualitative agree- 
ment in the case of temperature and very good overall 
agreement in the case of the Mach number distribu- 
tion. Figure 7 shows the comparison between the pre- 
dictions and experiment of the lower wall (hydrogen 
jet side) pressure. Ignition causes the pressure rise in 
the profile. Ignition seems to be delayed in the case 
of the predictions accompanied by a more pronounced 
pressure rise. 

High speed reacting flows such as the two cases stud- 
ied here are complex inspite of their simple geome- 
tries. The interactions between the different aspects 
of the flow such as turbulence, chemical kinetics, heat 
release etc. are very difficult to understand and, to 
a large extent, impossible to model accurately. Mod- 
ern day experimental facilities still cannot make com- 
plete measurements in such flows. Only mean values 
of temperature, velocity, pressure, species concentra- 
tions etc. are available, if any. Even then, the un- 
certainties associated with the data force one to ac- 
cept them only with certain reservations. Given that, 
there is almost never a chance for perfect agreement 
between predictions and experiment in all the areas. 
While the advances made in measurement techniques 
improve every day, the fruits of these advancements 
(ie. accurate measurements) are not realized imme- 
diately. As a result, today's computations will have 
only old data for validation (in the present case, the 
best data is already four years old) purposes which is 
certainly the case here. Unless more accurate experi- 
mental data with minimum uncertainties are available, 
the best a computational effort can hope for, in terms 
of validation, is probably what is seen here. 

CONCLUSIONS 

Computation of the flow fields of two high-speed, 
turbulent, reacting flow configurations involving finite- 
rate chemical kinetics for hydrogen-air combustion 
have been carried out. A two-equation ( k — e) turbu- 
lence model with compressibility corrections has been 
used. The predictions axe compared with available 
experimental data. Good qualitative agreement is 
present between computations and experiment. More 
detailed experimental data is necessary. 
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Table 1. H 2 — Air Reaction System 


No. 

Reaction 

1 

H + 0 2 — O + OH 

2 

OH + H 2 s* H 2 0 + H 

3 

0 + H 2 ^± OH + H 

4 

OH + O H-. -H 2 0 + 0 

5 

H + OH + M ^ H 2 0 + M 

6 

H + H -f- M ^ H 2 + M 

7 

O + O + M ^ 0 2 4* M 

8 

H + O + M ^ OH + M 

9 

H + 0 2 + M ~ H0 2 + M 

10 

OH + H0 2 ^ H 2 0 + 0 2 

11 

H + H0 2 ^ H 2 + 0 2 

12 

H + H0 2 OH + OH 

13 

O + H0 2 OH + 0 2 

14 

H0 2 + H0 2 5=^ H 2 0 2 4* 0 2 

15 

H + H 2 0 2 H 2 + H0 2 

16 

OH + H 2 0 2 H 2 0 + H0 2 

17 

H + H 2 0 2 ^ H 2 0 + OH 

18 

0 + H 2 0 2 ^ H0 2 + OH 

19 

OH + OH + M ^ H 2 0 2 + M 

20 

OH + OH v* H 2 + 0 2 


Species : H 2 , 0 2t H 2 0, OH, H, O , H0 2 , H 2 0 2 and 
N 2 (inert) 

M is a third body (all species included) 


Table 2. Conditions for coaxial jet experiment 



H 2 

Air 

Mach No. 

1.0 

1.97 

Temperature 

545 K 

1180 K 

Pressure 

0.112 MPa 

0.107 MPa 

Velocity 

1774 m/s 

1380 m/s 


1.0 

0.0 

fo a 

0.0 

0.246 


0.0 

0.545 


0.0 

0.209 


Fuel injector diameter=Q. 00236 m 
Lip thickness=0. 000725 m 
Nozzle diameter(air flow)=0. 01778 m 
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Table 3. Conditions for Burrows-Kurkov 
Experiment 



Hi 

Air 

Mach No. 

1.0 

2.44 

Temperature 

254 K 

1270 K 

Pressure 

0.1 MPa 

0.1 MPa 

Velocity 

1216 m/s 

1764 m/s 

f//* 

1.0 

0.0 

fo. 

0.0 

0.258 


0.0 

0.486 

f h 7 o 

0.0 

0.256 


Fuel injector height=0.004 m 
Duct height at inlet=0.0938 m 
Duct height at exit=0.1048 m 



U a j r = 1380 m/s 
= 1774 m/s 

2 

Fig- 1 


air 


1180 K 


T h = 545 K 

h 2 


Coaxial jet case 


P air = 107000 N/m 2 

P u = 112000 N/m 2 

H 2 

schematic 
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i-Kurkov expt. : schematic 
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ABSTRACT 


Computation of scramjet combustor flows is the 
main thrust of this work. Computations are done on 
flow configurations that resemble the reaction zone in 
the scramjet combustors. Compressible, reacting, tur- 
bulent flow solutions are obtained. A two equation 
(&-e) model with compressibility correction is used to 
calculate the flow field. A finite rate (9-species, 20- 
reaction steps) chemistry model for hydrogen-air com- 
bustion has been used. Computations are carried out 
using the Navier-Stokes solver TUFF. Predictions are 
compared with available experimental data and also 
those obtained by using the codes UPS and SPARK. 
A turbulence-chemistry interaction model is included. 
This model represents the first step in the effort to in- 
clude these interactions in high-speed reacting flows. 
A two-dimensional reacting mixing layer has been used 
for demonstration. 
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U velocity vector 

w n production rate of species n 

x streamwise coordinate 

Xj j th coordinate 

y transverse coordinate 

a, a* turbulence model coefficients 

ao,t*() turbulence model coefficients 

c Td , <7<jj turbulence model constants 

Sij Kronecker delta 

e turbulence energy dissipation rate 

e c compressible dissipation rate 

77 compressibility correction coefficient 

7 ratio of specific heats 

c j specific dissipation rate 

fi laminar viscosity 

fit turbulent viscosity 

1 / kinematic viscosity 

p density 

<Tk, <t € turbulence model constants 

Tij stress tensor 

Subscripts 

t turbulent quantity 


INTRODUCTION 

Propulsion systems for hypersonic applications are 
required to use the supersonic combustion ramjet 
(scramjet) due to the flight regime involved. Due to 
the advent of the National Aerospace Plane (NASP) 
proposed over a decade ago, there has been a consider- 
able amount of effort directed at developing a scram- 
jet based propulsion system for a number of years. 
Even though interest in the NASP is waning at present 
due to various reasons, the research effort in the area 
of high speed propulsion systems such a s the scram- 
jet still is going strong as evidenced by the enormous 
number of reports published in the literature and na- 
tional and international meetings. It is not within the 
scope of the present report to list and describe all the 
work going on in the area of high speed propulsion. 
However, references [1] to [11] are some examples of 
the vast amount of work in this area. The area of 
high speed propulsion, especially for hypersonic appli- 
cations, is a complex one encompassing a wide vari- 
ety of related topics such as inlet configuration, mix- 
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ing layers, mixing enhancement, combustor configu- 
ration, finite rate chemistry models and chemical ki- 
netics. Both hydrocarbons and hydrogen are consid- 
ered as fuels depending upon the application. For ex- 
ample, hypersonic waveriders using hydrocarbon fuels 
have been designed [1] for applications in the moderate 
hypersonic speed regimes. For Mach numbers of the 
order of 15 and above, hydrogen is generally considered 
to be the fuel of choice. The present work represents a 
computational effort in establishing a solution proce- 
dure for high speed propulsion in general and scramjet 
flows in paricular. 

The flow field in a scramjet is complex. It is 
turbulent and compressible involving high heat re- 
lease. The solution procedure should address all as- 
pects of the flow field adequately. It should be ca- 
pable of accurately modelling the turbulent field, tak- 
ing into account the effects of compressibility, address 
the changes associated with heat release and accuar- 
tely describe the chemical kinetics mechanism which 
is vital in high speed combustion. Also, the inter- 
actions between the distinct physical aspects of the 
flow such as the effect of heat release on turbulence, 
the interaction between turbulence and chemistry etc. 
must be properly addressed. Progress has been made 
in addressing these areas via accurate and realistic 
modelling in recent years [6]. However, a significant 
amount of work remains to be done before any real- 
istic claims for accurate models in these areas can be 
made. There have been remarkable advances made 
in the area of turbulence modelling, accounting for 
a variety of factors that affect the flow field. Com- 
pressibility correction models to account for the ef- 
fects of compressibility, near-wall turbulence models 
to deal with the transition from fully turbulent to zero 
turbulence, viscous dominated flow field near no-slip 
boundaries and modifications to models to account for 
flow curvature are some examples. A wide variety of 
turbulence models, ranging from the algebraic (zero- 
equation) models to Reynolds stress and large eddy 
simulation models are available (references [12] - [14] 
are some examples) depending upon the sophistica- 
tion and accuracy desired and the limits imposed by 
numerical solution procedures. 

Thermodynamic and chemical kinetic models [15, 
16] applicable to the scramjet flows have been under- 
going continuous improvements in recent years. Ac- 
curate modelling of thermodynamic variables as func- 
tions of temperature which are valid over a wide range 
of temperatures is an example. The design of the com- 
bustor is strongly dependent upon factors such as the 
mixing between fuel and oxidizer streams, presence of 
shocks in the flow field, boundary layer effects, flow 
separation, extent of chemical reaction within the com- 
bustor and so on. The numerical solution procedure 


should have the capability of addressing all of these 
factors while maintaining the required accuracy and 
robustness. Computational algorithms which are fast 
and accurate are being improved everyday and being 
adapted in present day numerical solvers applicable for 
a wide variety of flows including all speed regimes. In- 
spite of the existence of a wide variety of sophisticated 
and physically accurate models in many of the areas 
associated with high speed propulsion, it is not always 
possible to use the best models in a numerical sim- 
ulation due to the limitations imposed by computer 
resources, computational economy and adaptability to 
practical problems. Fast turnaround is a key factor 
in the success of any solver in practical problem solv- 
ing situations. Compromises must be struck between 
physical accuracy and computational feasibility. 

The present work represents the effort made to es- 
tablish a solution procedure for high speed propulsion 
related flows. The models chosen to represent the 
turbulent and chemistry fields reflect the compromise 
between physical accuracy and computational econ- 
omy mentioned above. The code chosen for the com- 
putations are the TUFF [17] code and the SPARK 
code [6, 9]. The solutions are compared with those ob- 
tained with the UPS [18, 19] code and available exper- 
imental data. The turbulence model chosen is the two- 
equation (fc — e) turbulence model with low Reynolds 
number modifications [12]. The compressibility effects 
are included via the compressibility correction model 
proposed by Zeman [20]. The fuel used is hydrogen and 
a 9-species, 20-reaction steps chemical kinetics model 
for hydrogen-air combustion [15] (Table 1) is used (in 
its abbreviated version). 

In the present task, major emphasis is on the turbu- 
lence model(s) used for reacting, compressible flows. 
Specifically, the interaction between the turbulence 
field and the chemical kinetics is of interest. The effect 
of turbulence on phenomena such as mixing between 
two streams, boundary layer etc. has been studied ex- 
tensively and it is an ongoing process. However, there 
is not enough work done in the area of turbulence- 
chemistry interactions in order to either emphasize or 
neglect the importance of such interactions. The main 
reason for this is the fact that it is extremely difficult to 
quantify these interactions experimentally. There have 
been efforts aimed at establishing pdf-based models for 
coupling the turbulence field with the thermochemi- 
cal field in numerical simulations. References [6], [21] 
and [22] are just a few of the many reports available 
in this area. There has also been a concerted effort to 
generate validation data bases through direct numeri- 
cal simulations which is still in its early stages. 

In the present work, on the thermochemistry front, 
the effects of temperature and species the effects of 
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temperature and species fluctuations are decoupled. 
For the temperature-turbulence interaction, a moment 
model is used [6]. The chemical species fluctations 
are accounted for using an assumed multivariate-/?- pdf 
model [6, 21, 22]. The thermochemistry models used 
here are independent of the model for the turbulent 
velocity field. These models require the solution of 
the evolution equations for ail the mean flow thermo- 
chemical variables, turbulence variables (it and either 
Lj or e), enthalpy variance and turbulent scalar energy 
(sum of all the species variances). The CFD code used 
is SPARK [9]. 

Mixing plays a major role in high speed combustor 
flows. The reaction zone is mainly confined to mixing 
layers that exist between fuel and oxidizer streams. 
Efficient mixing of fuel and oxidizer is very impor- 
tant for the design of combustor size. MLxing lay- 
ers ranging from subsonic to supersonic speeds have 
been studied extensively over the years. In the present 
work, a two-dimensional and axisymmetric compress- 
ible, reacting, mixing layer (hydrogen-air) is computed 
with the aforementioned codes. Experimental data 
from compressible, reacting mixing layers is still scarce 
which hinders the validation of the calculation pro- 
cedure. Three flow configurations are chosen for the 
study. The first two are chosen to demonstrate the 
solution procedure without the turbulence-chemistry 
interaction model and the last case is used exclu- 
sively to demonstrate the effect of this model. The 
first case is the well known Burrows-Kurkov experi- 
ment [23] in which hydrogen and vitiated air streams 
(two-dimensional) mix and react. The second case is 
that of an axisymmetric configuration [24, 25] where 
two coaxial jets (fuel and oxidizer) mix and react. The 
two-equation, (fc— e) turbulence model is used for these 
applications. Available data from the above two ex- 
periments are used to compare with the predictions. 
The last case considered is a two-dimensional, react- 
ing, mixing layer problem. Here, the two-equation, 
(fc—w) turbulence model [14] has been used. The gov- 
erning and secondary equations used in the computa- 
tions have all been described in detail in the references 
cited above. Only an abreviated equation set will be 
given in the present paper. The computations were 
performed on the supercomputers of NAS and NASA 
Ames Research Center. 


a tilde represents Favre (density-weighted) averaging 
and a simple overbar represents conventional time 
averaging. The averaged continuity and momentum 
equations are 


am 

dt dxi 
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dm dpU.Uj = dp dpu'-u'j dr^ 

dt + dxj dxi dxj + dxj { > 

where 




2 du k 

3 dx k ij 


( 3 ) 


with repeated indices indicating summation. 

In the two-equation turbulence model, the two tur- 
bulence variables are the turbulent kinetic energy ( k ) 
and either the dissipation rate e [12] or the specific dis- 
sipation rate u> [14]. The definitions of these variables 
can be found in the cited references. Boussinesq ap- 
proximation is used to obtain closure of the averaged 
equations. The Reynolds stress tensor is written as, 




- = Pi Sij - | pkSij 

dU\ dU) 2 dUl 

dxj + dxi 3 dxk 


( 4 ) 


where y.% is the turbulent/eddy viscosity defined as 



with C^=0. 09, or, with uj defined as, 


w = e/(Cs*) 

where 

c = ii+(£)! 

5 100 l + (^f) 4 

with R c =6 and 
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GOVERNING EQUATIONS 


fit is given as 


The equations used for computations are described 
in detail in references [6, 7, 10] and [26]. Only the 
forms of the modelled equations used in the present with 

study are given here. Density- weighted averaging is 
used to derive the mean flow equations from the in- 
stantaneous coservation equations. In the following, 
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is given below. 


«5 — Cs/3, Rk — 8, C/Q — “7n 


The modelled momentum equation, then, is, 

dpUi dpUiUj _ dp 2 (9p& 

+ dxj dx{ 3 dxj tJ 

- ^7 [(/* + *)$«] (12) 

The effects of compressibility are included via the 
model proposed by Zeman [20]. Here, the compress- 
ible dissipation terms are expressed as functions of the 
turbulent kinetic energy dissipation rate and the local 
turbulent Mach number. The compressibility effects 
are represented by a component of the dissipation rate 
(e c ) given as 

e c = K c e 
K c = rjF(M t ) 

2 k 

M t = ^ (13) 


where a is the local speed of sound and F(M t ) is a 
function of the local turbulent Mach number (M t ). 
F(M t ) is given by 


dpu dpuUj 

dt dxj 


- C 6 pw 2 + <7d—D k u 
k uj 


F(M t ) = 1 — exp[— ( 




— 0, M t < M io (14) 

with M* o =0.1 and 77=0.75. The modelled turbulent 
kinetic energy equation is [6, 12] 

dpk dpkUj _ , 

~dt + Ox, = p k~P ^( 1 + K c) 

d r/ _ fit . dk . 

+ ^ +f kW < 15 > 


where $ is either e or ujk and 

P k = 


— 77 — 77 dUi 

pUi ** ^7 


The modelled e-equation used (no compressibility 
corrections) in the present analysis [12] is given below. 

tr + : %7 i = IC ' Pt ~ C ' J, < l 

+ < 17 > 

where P t is the production term in the turbulent ki- 
netic energy equation. The modelled ^-equation [14] 



d . fl t dc 

(18) 

where 

dk du 

fcw ~ - 

(19) 

and 

1 <*o + 7^ 

a ~ 2a* l + ^r 

(20) 

The turbulence constants are ao=0.1, R^— 2.2, 

^=1.67, Ci = 1.44, C 2 = 1.92, <r fc = 1.0, <r £ = =1.3, 

Pr— 0.72, Pr t =1.0, Sc=0.22 and Sct=1.0. and 


< Td — 0.0, Dku < 0 

= 0.3, D ku > 0 



The mass-averaged total energy can be written in 
terms of the total enthalpy as 

E — H — \ (21) 

P 

The correlations between the fluctuating velocity and 
the scalar fluctuations cire modelled using a gradient- 
diffusion hypothesis. A typical model is of the form 

- (22) 

where is a coefficient which, normally, is a constant. 
For <f> = / n (n represents the species) , — Sc t} and 

for the static enthalpy, (<j> = h), cr# = Pr t . Using the 
above definition, and omitting the body force contri- 
bution, the time-averaged and modelled energy equa- 
tion [6] is 

dpE , dpEUj d _ z TTT!\~ 

~W + ~d^T = Yx 

, & u fi , Pt . dh ii, dk 

+ dxj ^ Pr + Pr t ^ dxj + ^ + cr* ^ dxj 

where <7* comes from the turbulent kinetic energy 
equation. The modeled species continuity equation is 

The modelled form of the mean species production 
rate due to chemical reaction (ti> n ) is given, for a finite- 
rate system involving L reaction steps and N species, 
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in the following general form: 
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Wn = M„ ]T(l/", - v' nl ) X 


(25) 
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exist. The factor m represents the effect of tempera- 
ture fluctuations on the mean reaction rate. The tem- 
perature variance is calculated from the enthalpy vari- 
ance [6]. 

For reacting flows involving multiple chemical 
species, a multivariate /?-pdf model is used to ac- 
count for the effects of the scalar fluctuations on the 
species production rates [6]. Relevant features of this 
model are briefly outlined below. The parameters of 
the multivariate /?-pdf for the iV-scalar mixing process 
(Ply •■*,/?#) are functions of the mean mass fractions 
f n and turbulent scalar energy Q : 


where, v‘ iX and are the number of molecules of the 
scalar s involved in the I-th reaction step in the forward 
and backward directions, respectively. The forward 
and backward rate-constants of the reaction / are given 
by kfi and k m respectively. 


kj, = A,r 4 ' exp[ ~^] (26) 

where A\ , 6r and Ta, are numerical constants specific 
to the given reaction step L k^i is determined from the 
equilibrium constant for the /- th reaction step and k/i. 

Turbulence- chemistry interaction model 


where 


Pn — fn( q 1) 

(30) 

N N 

s = E /« . Q = £ /»/«• 

n = l n = 1 

(31) 


Decoupling the effects of temperature and species 
concentrations, the mean species production rate can 
be written as, 


L 

= M n ~ V nl ) * 

/=! 


(32) 


Practical combustor flows involve multiple scalar 
mixing and reactions. In order to calculate the mean 
species production rate using the probability density 
function (pdf) approach, a model for the joint pdf of 
temperature and various species is required. Such a 
joint pdf is difficult to specify and as a preliminary 
step we effect the simplification that temperature and 
species fluctuations are uncorrelated. This permits us 
to deal with the temperature and species fluctuations 
separately. Equation (26) can be written as 

*„ = (27) 

Assuming that < 1, the term (T+T n ) can be ex- 
panded in a series and the resultant modified reaction 
rate term is written as, 

kj, - (1 + m) A,T b ' exp[— h-] (28) 


N 1 V 

5=1 5=1 

where 

5=1 5=1 

In the above, angular brackets represent conventional 
time averaging. The modelled expression for Ijj using 
the multivariate /?-pdf is [21], 

N v rl rni 

= n n^* + ^ ri( B + ( 33 ) 

5=1 r=l p = l 

Similarly, the expression for hi is, 

JV n, 

hi = n m + u "< - r )/ + n > - p) ( 34 ) 

5 = 1 P=1 p~ 1 


where 


- = [(fc-l>(| + ^) + (29) 


Terms of order higher than two in are neglected 
from the series expansion in the present analysis since 
accurate models for these higher order terms do not 


where 

* = 0i+ft + --- + Av (35) 

The production of turbulent scalar energy due to 
chemical reaction can be decomposed as 

N N N 

Wnfn = ^ ™nfn ~~ ) ] w n fn- (36) 
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Only the first term on the right hand side of the 
above equation needs further modeling and details 
can be found in reference [6]. The turbulence- 
thermochemistry interaction models require the en- 
thalpy variance and the turbulent scalar energy distri- 
butions. The modeled equations for enthalpy variance 
(h"h") and turbulent scalar energy(Q) are of a form 
similar to that of the turbulent kinetic energy: 


dpg 

dt 


dpgUj 


= - n 

Ptx dg 


dxj 
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(37) 


For g = h n h N ) G = h, ^=0, cr — Pr and <r g = Pr t . 

9 = Eli /£/». G 7 /«, *= 2Eli W" = 
■Sc and = Sc t . The dissipation term in the above 
equation is assumed to be 


€ 9 ~ ^ 9 — CgCiLjg 


(38) 


Solution of the modeled equations 


The equations are discretized and integrated in 
space and time to obtain steady state solutions us- 
ing either the finite-volume based numerical solver 
TUFF [17] (for the first two test cases) or the fi- 
nite difference solver SPARK [9] (for the last case). 
The TUFF code contains many desirable features for 
the computation of three-dimensional, hypersonic flow 
fields. It has non-equilibrium, equilibrium and perfect 
gas capabilities along with an incompressible option. 
It employs a finite-volume philosophy to ensure that 
the schemes are fully conservative. The upwind in- 
viscid fluxes are obtained by employing a new tem- 
poral Riemann solver that fully accounts for the gas 
model used. This property allows the flow field dis- 
continuities such as shocks and contact surfaces to be 
captured by the numerical scheme without smearing. 
Total Variation Diminishing (TVD) techniques are in- 
cluded to allow extension of the schemes to higher or- 
ders of accuracy without introducing spurious oscilla- 
tions. The schemes employ a strong coupling between 
the fluid dynamic and species conservation equations 
and are made fully implicit to eliminate the step-size 
restriction of explicit schemes. This is necessary since 
step-sizes in a viscous, chemically reacting calculation 
can be excessively small for an explicit scheme, and 
the resulting computer times prohibitively large. A 
fully conservative zonal scheme has been implemented 
to allow solutions of very complex problems. The 
schemes are made implicit by fully linearizing all of 
the fluxes and source terms and by employing a mod- 
ified Newton iteration to eliminate any linearization 
and approximate factorization errors that might oc- 
cur. Approximate factorization is then employed to 


avoid solving many enormous banded matrices. As 
mentioned before, the options for turbulence models 
include both zero and two equation models (both k — e 
and k - u/). For more details about the solution pro- 
cedure the reader is directed to the reference cited 
above [17]. The SPARK solver is a finite difference el- 
liptic solver using a hybrid McCormack scheme. More 
details about the solution procedure can be found in 
references [6] and [9]. 

RESULTS AND DISCUSSION 

The computational effort is divided into two parts. 
The first part deals with the demonstration and valida- 
tion process for the solution procedure by application 
to high speed turbulent reacting flows involving finite 
rate hydrogen-air chemistry. As mentioned before, the 
Navier Stokes solver TUFF has been used for the com- 
putations. Two reacting flow configurations have been 
chosen for this part. The first test case considered is 
the Burrows-Kurkov experiment [23]. The flow con- 
figuration is two-dimensional. A schematic diagram of 
the configuration is given in figure 1. No-slip walls 
bound both the upper and lower regions (y=0 and 
y—ymax)> The lower wall is inclined to form an expan- 
sion surface. Hydrogen is injected into a vitiated air 
stream. The two streams mix and react downstream of 
the injection location (inlet). The hydrogen stream is 
injected at a velocity of 1216 m/sec and a temperature 
of 254 K. The airstream comes in at a speed of 1764 
m/sec and a temperature of 1270 K. Full details about 
the flow parameters and geometry are given in Table 2. 
In this case, the reference length used in the hydrogen 
jet width at inlet, h (h=0.004 m). A constant turbu- 
lence intensity level is used for arriving at the initial 
distribution of turbulent kinetic energy and the dis- 
sipation rate. A 13-step, 8-species H o — Air reaction 
model (Table 1) has been used for the finite-rate chem- 
istry system considered here. The grid size is 81 X 121 
(81 grid points in the axial (x) direction and 121 grid 
points in the transverse direction). The total length 
of the solution domain is 0.356 m (x/h=89). Avail- 
able inlet data have been used for the first-plane pro- 
files which improved the predictions remarkably over 
the solutions obtained with uniform profiles. The so- 
lutions are compared with the available experimental 
data at this location (exit) in figures 2-4. The solu- 
tions carried out with the space marching PNS code 
UPS [19] using the Baldwin-Lomax turbulence model 
are also given for comparison. 

Figure 2 shows the comparison between the pre- 
dicted distributions of the species mole fractions and 
the corresponding experimental data. As seen in 
these figures, there is excellent agreement between the 
TUFF predictions and experiment. The predictions by 
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the UPS and SPARK codes do not agree very well but 
still there is very good qualitative agreement with the 
experimental data. Figure 3 compares the predicted 
profiles of exit plane total temperature with the ex- 
perimental data. There is good qualitative agreement 
here. Figure 4 shows the comparison between the pre- 
dictions and experiment of the lower wall (hydrogen jet 
side) pressure. Ignition causes the pressure rise in the 
profile. Ignition seems to be delayed in the case of the 
TUFF predictions accompanied by a more pronounced 
pressure rise. A point to note here is the fact that the 
SPARK, UPS and TUFF codes produced nearly iden- 
tical results for nonreacting mixing layers which is not 
the case in the present reacting flow simulations. 

The second case is that of coaxial jets [24, 25] where 
a hydrogen jet flows (inner jet) coaxially with an 
outer vitiated air (mass fractions: oxygen=0.246, wa- 
ter=0.209 and nitrogen=0.545) jet. A schematic of the 
flow problem is given in Figure 5. The two streams 
are, air (17=1380 m/sec, T=1I80 K with p=107000 
N/m2) and hydrogen ({7=1774 m/sec, T=545 K with 
p= 112000 N/m2)., The air stream is supersonic with a 
Mach number of 1.97 and the hydrogen stream Mach 
number is LOO. The inlet mean velocity is assumed to 
have a step profile with the two jets having uniform 
speeds at the specified values (no experimental data 
available). The velocity in the lip region of the inner 
jet tube wail (finite wall thickness) is assumed to be 
zero. The inlet temperature profile is derived based on 
the experimental data given for a location just down- 
stream of it (shown later). The inlet species mass frac- 
tion distributions are also chosen based on the exper- 
imental data provided at the same downstream loca- 
tion. The models for turbulence and chemistry are 
identical to the ones used for the first test case. A 
81 X 91 grid (81 points in flow direction, 91 points 
in the radial direction) was used for the calculations. 
The inner jet/tube diameter (77=0.00236 m) is used 
as a reference length. The total length of the flow do- 
main is equal to 43.1 D. The outer boundary (radial) 
of the flow domain is taken to be at y=17 D. A more 
detailed description of the flow parameters is given in 
Table 3. The region outside the limits of the air jet is 
assumed to be still air at a temperature of 273 K. The 
two-equation (£ — e) turbulence model is used along 
with the finite rate H^-Air chemistry model mentioned 
above. In all the figures shown for this case, y refers 
to the radial distance measured from the axis of the 
coaxial system of jets. 

Figures 6-7 show the results of the computations. 
Figure 6 shows the computed and experimental distri- 
butions of species mole fractions. The figure is de- 
signed in a two-column format. The left side col- 
umn represents the inlet (first x-location) data and 
the right side column is the data at the exit plane 


(*/25=43.1 D). As seen in these figures, the inlet 
data agreement between the computations and exper- 
iment is not perfect, especially around the jet edges, 
and this might affect the computed distributions at 
downstream locations. The comparison between pre- 
dictions and experiment at the downstream location 
{x/D— 43.1 D) is good given the above mismatch be- 
tween the computational and experimental data at the 
inlet. The development of the reaction zone after igni- 
tion is not predicted well by the TUFF code whereas 
the SPARK code fares better. The experimental data 
indicates that the reaction zone (depicted by the water 
mole fraction distribution) spreads more quickly than 
the predictions indicate. There seems to be a discrep- 
ancy between the locations of peak reaction activity 
between the predictions (off the axis) and the exper- 
iment (closer to axis). However, there is very good 
qualitative agreement berween the data with the peak 
values of the reaction products predicted very well. 
The flow domain was seen to have a wave-like structure 
as shown by the predicted profiles. Figure 7 shows the 
comparison of static temperature data. The agreement 
between predictions and experiment is good qualita- 
tively displaying similar trends. The uncertainty as- 
sociated with the accuracy of the experimental data 
is unknown. There are considerable differences be- 
tween the data presented by the two references [24, 25], 
especially in the temperature profiles. Overall, there 
is good qualitative agreement between the predictions 
and experiment. 

In the second part, the turbulence-chemistry inter- 
action model is demonstrated using a two-dimensional 
mixing layer configuration. A schematic of the flow 
problem is given in figure 8. The two streams 
are, air ( U =1606 m/sec, T =1600 K with 
fit 2 —0.0, fo 1 = 0.267 and /yv 3 =0.733) and hydrogen 
( U =1250 m/sec, T =254 K with //f a = 1.0, /o 3 = 0.0 
and //v 3 = 0.0). The two streams are supersonic with 
the air stream Mach number of 2.07 and hydrogen 
stream Mach number of 1.03. The inlet mean veloc- 
ity is assumed to have a hyperbolic tangent profile. A 
constant turbulence intensity level is used in the free 
stream for arriving at the initial distribution of tur- 
bulent kinetic energy and the specific dissipation rate. 
The pressures are matched between the two streams 
(P =1 atm.). A 13-step, 8-species H 2 — Air reaction 
model (Table 1) has been used for the finite-rate chem- 
istry system considered here. A 101 X 81 grid (101 
points in flow direction, 81 points in the transverse di- 
rection) was used for the calculations. The length of 
the flow domain is 0.25 m and its width is 0.05 m. In 
the figures the title g refers to calculations that include 
the interaction model. 

The difference made by the temperature fluctuations 
model is shown in figure 9. Here the factor (1 -h m) is 



plotted as a function of y for the first three reaction 
steps at the exit plane of the flow domain. The need 
for including the temperature fluctuations is evident 
in the figure. The forward reaction rate changes by 
as much as 400 percent (reaction step 1) due to the 
efFect of this model. The species production rate is a 
strong function of the reaction rate ( 29) and hence 
will be affected by the inclusion of this model. The 
effect of the turbulence-chemistry interactions model 
on the species production rate is shown in the next few 
figures. Figure 10 shows the production rates of the 
major species involved (H 2 , 02 1 ^ 2 0, OH) as func- 
tions of y at the exit plane. Computational results 
obtained with (k-cj-g) and without (k-cj) the effect of 
these interactions are shown in the figure. It is seen 
that the interaction between turbulence and chemistry 
has a significant efFect on the net production rate of 
chemical species. While it is not conclusive whether 
the turbulence-chemistry interactions increase or de- 
crease the production rates of individual species, it 
is important to note that these interactions do have 
a significant efFect on progress of the chemical reac- 
tions. Figure 11 shows the efFect of the interactions on 
the species distributions. The efFect seems to be more 
pronounced in the case of the products (H^O, OH) 
than the primary reactants (# 2 , 02). The changes 
are larger in regions of higher gradients (on the fuel 
stream side for H 2 and air stream side for 0o, for ex- 
ample). 

This effort represents a significant push in the right 
direction in this very important area of turbulent re- 
acting flows. The model is simple to use and is easily 
adaptable to other types of turbulence closures such 
as the Reynolds stress models. More improvements 
can certainly be done in terms of using more realistic 
models for the temperature and concentration fluctua- 
tions such as a pdf-based model for temperature effects 
(as opposed to the moment model employed here), a 
joint-pdf model to represent the effects of both the 
temperature and the concentration fields etc. How- 
ever, in order for the model to be useful for practical 
applications, it must be kept in a form which promotes 
ease of use as well as computational economy. This is 
the main reason for resorting to a two-equation level 
turbulence model in the present work. A concerted 
validation effort is crucial to ensure the success of the 
modelling efforts. One debilitating factor in this re- 
gard is the nonexistence of useful experimental data. 
It is extremely difficult, if not almost impossible, to ob- 
tain such experimental data using present day equip- 
ment. An alternative seems to be direct numerical 
simulations (DNS) which has shown promising signs 
in relatively simpler flow configurations. 

CONCLUSIONS 


A numerical solution procedure applicable to high 
speed reacting flows has been demonstrated. In the 
first part, computation of the flow fields of two high- 
speed, turbulent, reacting flow configurations involv- 
ing finite-rate chemical kinetics for hydrogen-air com- 
bustion have been carried out. A two-equation ( k — e) 
turbulence model with compressibility corrections has 
been used. The predictions are compared with avail- 
able experimental data. Good qualitative agreement is 
present between computations and experiment. More 
detailed experimental data is necessary. In the sec- 
ond part, a turbulence-chemistry interaction model is 
introduced and its effect is demonstrated via a two- 
dimensional turbulent reacting flow problem. The 
interactions seem to have a significant effect on the 
species production rate. More detailed analysis is nec- 
essary. 
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Table 1. Hn - Air Reaction System 


No. 

Reaction 

1 

H + O 2 -O + OH 

2 

OH + H 2 H 2 0 + H 

3 

O + H 2 - OH + H 

4 

OH + OH — H 2 0 + 0 

5 

H + OH + M H 2 0 + M 

6 

H + H + M ^H 2 + M 

7 

O + O + M 0 2 + M 

8 

H + O + M^OH + M 

9 

H + O 2 + M ** HO 2 + M 

10 

OH + HO 2 — H 2 0 + 0 2 

11 

H + H0 2 H 2 + O 2 

12 

H + HO 2 OH + OH 

13 

O + HO 2 w* OH + O 2 

14 

HO 2 + H0 2 H 2 0 2 + 0 2 

15 

H + H 2 O 2 ^H 2 + HO 2 

16 

OH + H 2 O 2 w* H 2 0 + HO 2 

17 

H + H 2 O 2 — H 2 0 + OH 

18 

O + H 2 O 2 — H0 2 + OH 

19 

OH + OH + M H 2 O 2 + M 

20 

0% + OH u* H 2 + O 2 


Table 3. Conditions for coaxial jet experiment 



a 2 

Air 

Mach No. 

1.0 

1.97 

Temperature 

545 K 

1180 K 

Pressure 

0.112 MPa 

0.107 MPa 

Velocity 

1774 m/s 

1380 m/s 


1.0 

0.0 

fo 3 

0.0 

0.246 


0.0 

0.545 

0 

c* 

0.0 

0.209 


Fuel injector diameter=0. 00236 m 

Lip thickness=0. 000725 m 

Nozzle diameter(air flow)=0. 01778 m 


Species : #2, 02> H 2 O, OH , H t 0, HO 2, H 2 O 2 and 
N 2 (inert) 

M is a third body (all species included) 


Table 2. Conditions for Burrows-Kurkov 
Experiment 



h 2 

Air 

Mach No. 

1.0 

2.44 

Temperature 

254 K 

1270 K 

Pressure 

0.1 MPa 

0.1 MPa 

Velocity 

1216 m/s 

1764 m/s 


1.0 

0.0 

fo a 

0.0 

0.258 


0.0 

0.486 

Of? 0 

0.0 

0.256 


Air (vitiated) 


U • = 1764 m/s 
air 

T air = 1270 K 



Urr = 1216 m/s 
z 

To = 254 K 
z 


Fig.l Burrows-Kurkov experiment : schematic 


Fuel injector height=0.004 m 
Duct height at inlet=0.0938 m 
Duct height at exit=0.1048 m 
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Fig. 6 Coaxial jet case : species mole fractions 








Fig. 8 Mixing layer (2-D) - schematic 
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Introduction 

O NE of the major areas of current interest where compu- 
tational fluid dynamics (CFD) is used extensively is the 
development of advanced air-breathing propulsion systems for 
hypersonic vehicles. A hydrogen- fueled supersonic combus- 
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comparison between the predictions and representative experi- 
mental data. 6,7 Here, C 6 is defined as 

= dd U 1 +JJ 2 
s dx Ui - U 2 

and Q 0 is its value for incompressible flow. KEPS1 refers to 
the predictions with compressibility correction included, and 
KEPS2 refers to the predictions without the correction. RS 
refers to the Reynolds stress closure predictions. 9 The scatter 
among the experimental data is indicative of the nature of 
available data. The effect of compressibility on the mixing- 
layer growth rate is to reduce the growth rate with increasing 
convective Mach number. The two-equation model without 
the compressibility correction does not predict the reduced 
growth rate with increasing compressibility. However, the 
compressible two-equation model does predict this trend very 
well. The close comparison between the Reynolds stress clo- 
sure and the k-e model is important because for complicated 
flowfields, such as that in the scramjet combustor, the higher- 
order Reynolds stress closure may be very expensive to use. 
The predictions agree well with experiment within the scatter in 
the available data. It must be pointed out that there are other 
compressibility correction models, similar to the one used 
here, available today. However, a detailed comparison be- 
tween such models is beyond the scope of this report. 

Mixing between fuel and air and the ensuing chemical reac- 
tion is the main focus of study for a configuration such as the 
scramjet combustor. The presence of turbulence and its effect 
on the flowfield, especially the mixing aspect of it, is an impor- 
tant part of such studies. A representative mixing dominated 
reacting flow (hydrogen and air) is considered here. A seven- 
step, seven-species H 2 -Air reaction model (Table 1) has been 
used for the finite rate chemistry system considered here. 

The test case studied is the reacting coaxial jet problem 
(Figs. 3) for which the experimental data was obtained by 




b) Centerline Hi distribution 

Fig. 3 Reacting flow: coaxial jets. 


Evans et al. 8 In the experiment, hydrogen at a temperature of 
251 K was injected (Mach number = 2.0 and velocity = 2418 
m/s) along the axis of a supersonic jet (Mach number 1.9, 
velocity = 1510 m/s) of vitiated air (temperature = 1495 K). 
The species mass fractions in the airstream are nitrogen = 
0.478, oxygen = 0.241, and water = 0.281. The fuel injector has 
an inner diameter of 0.6525 cm and outer diameter of 0.9525 
cm. The outer jet (air) diameter is 6.53 cm. The predicted 
centerline distribution of hydrogen and the profiles of the 
major species mass fractions (H 2 , 0 2 , H 2 0, N 2 ) at three axial 
locations are compared with experimental data in Figs. 3. In 
these figures, D is the fuel injector outer diameter. The predic- 
tions agree well with experiment in the region where the mixing 
effects dominate. The discrepancy between the predictions and 
experiment immediately downstream of the injector exit may 
be due to the fact that the conditions set at the inlet location 
(boundary condition for computations) in the calculations 
may not match the exact experimental conditions. These con- 
ditions are unknown and, hence, could not be used for the 
calculations. The predicted profiles of water vapor, which is an 
indication of the extent of reaction, agree reasonably well with 
the experiment. The location of the peak in the profile is 
farther into the airstream for the predictions than the experi- 
ment indicates. This may be because the interaction between 
turbulence and chemistry is not fully accounted for in the 
calculations. As mentioned earlier, the calculations account 
for only the effect of temperature fluctuations on the reaction 
rate. In addition, correlations of order higher than 2 were 
dropped from the model. Also, the experimental data itself 
may not be very accurate. Given that, the compressible turbu- 
lence model used here performs very well in predicting the 
turbulent reacting flow. 


Conclusions 

A two-equation turbulence model ( k-e ) has been modified to 
be applicable for compressible flows by adapting a compress- 
ibility correction model. Computation of compressible mixing 
layers indicate that the decrease in the growth rate of the 
mixing layer with increasing convective Mach number is well 
predicted by the model. Comparisons of the predictions agree 
very well with available experimental data and the predictions 
of a compressible Reynolds stress closure. Preliminary studies 
of reacting mixing problems involving dissimilar gases indicate 
that the model is well suited for application to such flows. The 
computations use a model for the effect of temperature fluctu- 
ations on the reaction rate in the finite rate chemistry model. 
Further improvement to the model, especially the effect of 
turbulent fluctuations in species concentrations on the produc- 
tion rate, is necessary. 
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